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PREFACE 


The present preprint is a continuation of the description, begun 
in (^), of the deicriptlon of the . libi-ary of subprograms^ arising in 
the process of creation of a packet of applied programs for the cal- 
culation of navigational parameters of artificial tteth ■etellltee of the 
Earth (^,8 ), wherein are contained the subprogx'ams having indexes 
from F to I, In the subprogram**; for the calculation of the right 
members of the iW8 enuatinns of motion, presented in Chapter I, 
the subprograms described in (5) are utilized. Here there operate 
also the principles of scaling the dimensional quantities, determined 
Ui (I). The constants and scale iviciors for the subprograms of chap- 
ters 1, 3» and k are routed ^ to the domain COMMON by conversion to 
the subprogram CONCTi. The system of coordinates determined in (3) is 
employed. 

The suoprogram for Integration of the system of differential 
equations by the method c ? Adams, described in Chapter 2, on the 
other hand, is sufficiently autonomous and may be used for integratiun 
of any system of ordinary differential equations. 

In Chapter 3 are presented subprograms ensuring the fixation 
of the attainment in the process of integration of the assigned 
values of different functions from the solution of systems of 
differential equations; the minimum and maximum of an arbitrary 
continuous function from the solution as a function of an independent 
variable; the exit of the AES at the ascending node of the orbit; 
the minimal and maximal altitudes of the' AES above the surface of 
the terrestrial ellj.psoid. 

Chapter k contains subprograms for the computation of the values 


iii 
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Of various functions from th® parameters of motion of the AES , 

The subprogram ROOT4(Ill), used for the computation of the 
moments of entry of the AIS into the umbra of the earth, is 
intended for the calculation of the roots of the algebraic 
equations of the fourth, third and second degree, has a 
significantly independent character and may be used for the 
solution of other problems# 

The subprogram FA QRAV (F 03) was written by E, E. Ryazanova. 

For the computation of the geomagnetic parameters B, L the subprograms 
BL, INVAR, LINES. STAR, CARMEL, INTEG, NEWKAG, ASIN (I 04) are 
utilized, submitted through the courtesy of Yu. N. Gal'perin and 
V, M, Sinitsyn, The author of the remaining subprograms is the au- 
thor olf the px^eprint. 

The author wishes to express his thanks to E. A, Chistyakova 
for assistance in the editing of the texts of the subprograms for 
publication, to L, V, Saytseva and V. F. Smirnova for help in the 
preparation of the manuscript. 
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^UBFi?0GHAMi5 FOR INTSaHHTl^«G THE EnUATlCiKS 
CF MOTION OF GATSLLITEO. FORTRAN FOUR. 

V. I. FrokHorenko 

Institute of Space Reosiirch, Academy of Sclcncss USSR, Moscow 


Chapter I. 


Right membaro of the system of ©euations of motion o<' 
^ (index F) 


I.l. Enuatlon of motion of the A18 . Various models. 

We will writ© the differential equations of motion of the 
AES in an absolute system of coordinates in the general case in 
the following form; 

x-]^,y.Vy,z.v,, 

^ ( 1 . 1 ) 

Vy t Vy f 

m 


In the nroenwlch relative system of coordinates these equations 
have tho form 

»,« o 4 jc *2ia^Vy f 

%*^y ^^sPx * * ^L^y » 

whoro tho projectiaiun of acceleration, determined by the influence of 
the con esponding forces, are: 


•Numbers in the margin indicate pagination in the foreign text," 


A/vftf * -normal gravitational field of the earth 

resistance of the earth's atmosphere 
-gravitational anomalies 

-gravitational perturbations of the moon, and sun 
Aji — pres.mro of light 

o)Jx,a»*sy — ^^centrifugal force 

2<OjWy,2Wjt», — force of Coriolis 

0>^ —angular velocity of the earth's rotation 


For the calculat I'n of the right mombcrB of the system of 
equations of moti.jn of AS8 one has the collection of subprograms: /G 

FNQRAV, FAGRAV, FATM, FGHh/., FUGHT, each of which allows for its 
auinncmont in the right membare of tre equations of motion of the 

AfiS*. 

Theue oubprograms have a subsidiary nature; from them it is 
possible to construct the subprogram for ..alculatlun of the right- 
hand sides for *r system of equations, with this or that degree of 
completion for the described motion of the In (4) were intro- 

duced the indexes KG, KG, KA, KG, KL and table 2,1, permitting the 
regulation of the variants of the system of forces, acting on the 
AES (and of the system of coordinates, in which the motion of the 
AES is analyzed), 

Table 1,1 is a repetition of table 2,1 from work (4), 


We recall also the dosignatii n of inuexes. 

The index KS characterizes the system of coordinates (pns ible 
values 1,2), 

To each of the forces acting is attached its index; 

KG — force of the earth's attraction (possible values 0,1,2,3»4); 

KA — atmospheric resistance (possible values 0,1, 2, 3, 4); 

KS — gravitational perturbation by the moon and sun Cpossible 
values 0,1, 2,3); 

KL — pressure of light (possible values 0,1,2). 

Giving to each of the above enumerated indexes the value de- 
termined, we give the determined model of forces, characterized by 
the five-valued index, consisting of tne values of the indexes 
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KG, IG, KA, KS, KL. 


^here it a tubprofraffl FOiOI, rtalitiBg all postibli variants 
of tbt aodilt of fore®, spaoifitd by tabl® 1.1. Howsvsr for aonerst# 
«o4tls it it possible to recommend to the user the creation of 
"truncated” subprograms, which may be obtained from the subprogram 
FORCI by means of discarding oonvarsions to those subprograms, 
which are not utilized in the given model. 


Treble 1.1 System of forces, considered in the equations of 
motion of the MB . Indexes KC, KQ, KA, KS, KL 



Systems 
p’f Coord. 

Grev, 

Field 

Atmo- 

•phare 

Grev. 

Perturbs^. 

tlon 

Presaure of 
Ll|ht 


KC 

KO 

KA 

KS 

KL 


1 

2 

3 

4 

5 

0 

oscular 
tlog elta, 

central 

not eofi- 
aldered 

not eon- 
■Idered 

not considered 

1 

GreenliAch 

r#lativs 

normal 

wlthoiit 

aalc*yar. 

from 

moon 

without calculatlcn 
of the earth's unbra 

2 

absolute 

sonal 

harmo- 

nica 

with 
calc. of 
long 

period var. 

from 

eus 

with calculation of 
the earth's umbra 

3 


iH^allei 

with calc, 
of long 
and short 
period var* 

from moon 
and aun 


4 


harmo- 

nlce 

22,30,40 

CIRA-72 




ae 




il 


\ 


The subprogram FG1100(FO6) ,may serve as an example of such a /8 
subprogram^ in this subprogram is realized the model of the motion 
of the ASS, corresponding to the following values of the indexes; 

KC equals 1 or 2, KC equals 0 or 1, KA equals 1, K3 ecuals 0, KL 
equals 0. 

In order that the designations of concrete subprograms-, of 
the right-hand sides, created by the user, reflect the model, 
realized in them, of the motion of the ABS, it is suggested that 
one construct identifiers of these subprograms according to this 
sort of principle; the first letter — F, after it a five- valued 
index, corresponding to the chosen variant of model from table 
1*1 (as this was don© in the case of FC1100). 
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1*? Culculation of thi norwal gravitati ml fiald of ths tarth and 
(for the case of the Greenwich systea of coordinates) translational 
and corioli® accelerations (FOl-FNGRAV) , 

‘ arc deteraiui^d; 

they art dependent on the influence of the normal gravitatl nal field 
of the earth in the right members of the equations of motion of the 
^ (1.1) in the absolute system of coordinates. If the motion of 
the A18 is calcuiuted in Greenwich system of coordinates (system of 
equations 1,2), then these sumso are determined; 


2, Structure, Thu subprogram FiJGRAV, General units; /CA00/.» , 
/CK/./^, /CAE/^, /CAEL/;j_, /COraEP/^. 


3. Conversion; CALL FNQHAV (KG, KG, Y, HC, RC, F) 

4, Initial data; KC — index, characterizing the system of co- 
ordinates (KC = 1 for Greenwich system of coordinates, KCsF for 
absolute evatam of coordinates;; 

KG— index, characterizing the gravitational field of the earth 
(for KG equal to zero there is considered only the central field) 

■*<%— main part, containing X, Y, ?, or x, y, z, 

V^, Vy, V^. 

3. Results; ££ 

HC — elevation of ASf above the surface of the earth; 

RC — modulus of the radius-vector of the AES; 

— main part of the ri^ht members of the equation of motiun of the 

AI8 , containing X, Y, S, V.r,Vv,V„ or A, ft Ktf. 

A X 4/ A y z 

6. Utilization of the groups COMMON; The constants of the 
groups are utilized; /CAOC| , /CR2^ , /CAE/ 2 , /°AEL{, /COMSP^ 

(see in (( 5)) . Nos. 2, 10, 11, 12, 15 table 2.1) 

7, Algorithm; 

a) absolute system of coordinates; 


liesignetl m. The term® 




k 


X'V,; . Y*Vy, Z-Yg, 

Vy ‘A.jyVy —oi„RX/r*-C(/>-l)llX/r* 
Vy •AjyVy — a.ty^Y/r*-C0-i)«r/r^ , 
Vg-AjyVji‘ -oi^Z/r*~C(D-S)KZ/r*, 


whar® Cs5/? «^Q(R/r)^, D«SZ^/r' , r«(X^+y^+''.^)^'^^* R— »v«rage radius 
or thf o«rth, «lphaQQ, alpha^p — paramoters of iht normal gravita- 
tional riaia of the earth; for alphaggsO there is coneidertd only 


the influence of the central gravitational field of the earth; 


b) Greenwich system of coordinates *, 


V 


X > 




V 


2 


y t 


Z mV. 


X » 


X 


“s* +2W,W^ 

t5y . A ;y + <0^ y - 2 W , t>. 


c) elevation h^ of the aes above the earth's surface ellipsoid 
is computed according to the formula h^sr-Cag-alpha^gVr^) , where 


a^, alpha— Beniimaj or axis and constriction of the general terrestrial 
ellipsoid. 


8, Text 


n° 
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$UStoi|T*Mi P»MlRAVtKC*Kt» 

fflMKNSirOM 

COMMON’ /cAOO/AOOtAfO 
COMMON /CAi/Af*A|. 
COMMON/OAfl./Ail. 
POMMON/f!l|Z/RZ 
»eoMMON <<eOMZf/OMZt(OMZK 

Rt>M«va)»V(|)4yCt)Ay(e> 

WaM/Nt‘ * 

NaSiRtilt) 

MC«t-AE*Atl.»l^ * 

Mi«NZ/N 
BaWl/Rt 
A*B«A0a 
BC*»0 * 


y > He * M » r > 


1 BCf *JI, #ARO*N4*i(fi#j0 
A*A*»ef</8.0f^W«,9Jf 

2 00 B .(vifO 

r«jAssB»-A«4MaJ 
9 p<«l)ay( 

PCOlKf f A)^Beti»NifAI 
OOTO <A»9>»|C|i 

r(9)!*rf9)A0MZKa)lf2$.p|iZf»y<A) 

!» RETtiNN' 

ENO 


i,3 €alcui.itlan of the Influsnc® of atmospheric resistance 
(FOr-FA'm), 

• • ♦ 

1, Deol(,mitian. The tar.iis *A4^ » AaH' ? Al^ • 

determioed, doi findlni^ un ^the influence^of the atiuosphci Ic resis- 
tance In the equMtiuns uf (l.l) or (l.T), 

r. structure. Subprot*ram FATK, 


General groups i /DSB/j., /COMS/^. 

5. Gonversion: CALL FATM (KG, Y, 1, F). 

If, Initial data*. KG— index, characterizing the system of co- 
ordinates 5 Y^— main part, containing X, Y, A, V^, Vy, Vj, or 

X, y, z, Vj^, V ,v^‘, r— density of the atmosphere; F^--fflain part of 

right me.ahers of aouatiuns of motion ofAlS X, Y, H, Vj^, Vy, 

or X, y, z, v„, v , v_ with the computati;*n of the influence of 

A J w# 

atmospheric resistance. 


6 , Use of the domain COMMON. Before conversion to the sub- 
program FATK in the group CCMMCN/BSB/SB it la necessary to addresr; 
the value of the ballistic coefficient in the system of units, 
obtained from the calculations, for this it is sufficient th«t the 
value of the ballistic coefficient, given originally in the system 
of units kg, m, s, be divided by the scale factor from the group 
COKMCK/CSSB/ESB (see in (5) No. IZ table 2.3). 
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from th# Kfoup C0KMCK/CCH::/CH'*, thtr« In utlUxtd thft con«tant 
(•if In (v) No, 14 t-Tbl# r.J) 

7* Algorithm 


Vr *AaV 


A’X • 


A/,Vx •-3B-PV 

»-aa‘PV 

Af^Vg *‘SB‘Py 


VynFy 


’ Aim., * 

m m 

/ ^ 
Vg •fg^A^wg 


<Tm 


where ^B«-balliatlc: co«f;ioi#nt 3B*0^,*Ft./Za 

X r* 


•di.itenuUmlffcti^ cof fficlant of air rRoiaiamt 


square oi too aid section, m— aaso of AM , 


I --nir aoimity 


yflan, y|!«h, yO.rm of the 

to the rtlr 

In theGre«-nwich system of coordinates 
In the absolute syotem of coordinates 


vector of tilaht sfeed relative 

,H 


:n¥7\\H 




•»V, +K 


j: 


vj^. , V^. v^, 

vr"‘>i^n»rf»{>-«b*f, >c-vz- 


F.., Fv. F-. 
orU.F). 


«-ot!i»-r tormn iuthr ri^ht iuo-.ibcrn of eouHtlons (l.l) 


o, fort 


i 

3 


sm^ooutine rATMfKC»y.^«r) 
COMMON /BCB/SB 
COMMON /«OMZ/OMZ 
dimension yc^) » v<»> 
O'O i J*i*3 
V(J)«V( 

80TO(Z»3) 'EC 

V(i>«vfi)»CNz«y(2» 


v(f )«v(B>'*OMr*vM> 

t HsO 
00 A 

4 WaMfyC 

N«SB*Nef09TCM) 

DO * 

betorn 

. .two. 


i./». Calculation of anoanlies of crnvlt* t^onal field of the earth 
(F 05 -FtGn'.¥) 

1. i)asignatl 'n. The terms delts„J/^, deltsQVy, deltaj,^ in 

the right membera of the system of enuatims (1.1) or (1.?), derived 
from the influence of anomaliea of the earth’s gravitational fiild. 
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structure. Subprogram FAGHAV, /Ig 

Thtr® are utilized the exterior subprograms; DSGF, DEG3, DKG4 (EDI) 
General group; /HAD/^. 

3. Gonveroion; ChLL FiiGHAV (Y, XG, P, KG, NM). 
l|.. Initial data. 

The main part Y^, containing X, Y, 'A, Vyj Vy, ^r x, y z 

Vjj, Vy , Vg*, the main part XG^, containing x, y, z; 

— main part of ripht merabora of equatiuna (1,1) or (l.T), 

containlnc the termsi determined by the influence of other factors; 

KG— index, determining the nature of the considered anomalies of the 
gravilntii.nal field of the earth (possible values; ?, 3, 4); 

For KGss? one considers only the zonal harmonics in the breakdown 
of the gravitational field of the earth, for KGs5 there are 
considered the zonal, teaseral, and sectorial harmonics, for KQ~h 
one ctmsidorH only the harmonics FF, 30, ifO; 

NM-~numb 0 r of considered harmonics (KK less than FF), 

Note 1, For KG»r or 3 it is necessary first to turn to the sub- 
program CONGR(AOF) for the addresMes In the group CCKMCN/BC0KQR4 

coofficients of the gravitational field of the earth. For KC=A 
one uses the coefficients ^rom the block COMKGN/GAFF/if (see (3) 

No. A table 2,1), value NM is not used here. 


Note F, The main part XG is used only for KQ=3 or 

• • • • 

3. Results: Main part F^- containing new values X, Y, X, Vv» 

• • ,,,,,,0 A 

Vy, V,. or X, y, z, v^^, Vy, v^ with calculation of the influence of 
the anomalies of the earth's gravitational field, 

6. Use of the domain COMMON, In the block COInMON/RAD/RC,Rl 
one nddres..ee the values RC=r and Rl=r^. 

7. Algorithm; 


; =-Ag,A/r-Ag„Zr/>rr, -Ag^y/r, 
A<,lV=-Ag^r>-Ag„Z<?rr, +Ag^Vr, 
Atf *-Agr^+Ag„ r,/r . 


,e’ 


■> 
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V^rFjj+delta^V^, VY=Fy+dsltagVy, V.,KF.,+d9ltiiQV.;,, wher® r=(X^+Y^+2^)^, 

p p 1 '* 

(X +y )^, Fy» P^^-cofflpontnfcs of othor terms Ln th^.' right 
membsro of equations (X.l) or (1.F)| 

delta gj,-- 'radial component of vector delta g, acceleration due to 

the Influence of anomaliaa of the earth’s gravitational field, 
delta gjjj— meridional component of the vector delta g, 

delta gj^--proJecti n nf the vector delta g on the normal to the 

plane of the meridian. 


Note, The components of the vector delta g, for some or some 
other recommended anomalies of the earth’s gravitational field, 
respectively, are calculated with the subprograms w th index SOI; 
i/SG ?, DSG 3 or i#SG 4 see in (3) p, 6.1 (there is introduced also 
the algorithm for the computation of the components of vector delta g) 


8, Text, 


1 

2 


3 

4. 


SU&RCUTINE P46P*VfX,K6»r,JV,K^> 
cimension y ( 6 ) » r * 6 si I C6C 3 > 
CIMENS!C*« XC(3> 

COMi^ON/RAC/Rl . R 

y2*xil)*xa) 

V.t*X{2)*X<2) 

22*X(3»*X( 3) 

Ri23X2*v2 
R1*S3RT (’*12 > - 
R2sel2*Z2 
RsSsR?<R2i 


NV* J V-1 
V^sR«Rl 

T S ( 1 » 2 • 3 ) , ^ V 
CAtt CE52( X ,06 ikM) 

GO 70 *> 

CAtL CESO’ XG.DG.nr) 

0 C TO A 

CALL 0EC*‘'^Oi0G) 
r(5)sr(5)-<CCU)*X(2)>/R- 


{06(2UX(1)*XC3) )/W-(0C(3>*X(2} )/R1 
(0G(2)*X(2}«*X(3) 3>*X^ 2 )>/Rl 


r(6)xF<6>-<CCCl)»X(3)>/R*l06<2J#Rl>/R 


RETURN 



l.l Calculation of the gravitational perturhatlons connected 
with the influence of the moon or sun (F04-FGRSS), 

• • 9 

1 , Designation. The terms delta^V^, delta^Vyi delta^V^ 

in the right members of the equations of motion (1.1) or (1.2), 
dependent on the gravitational perturbation of the moon or sun. 

2 . Structure. Subprogram FGRSS. 

5. Conversion: CALL FGRSS (GDR 3 , Y, XS, RS, F). 

h. Initial data: GLR 3 = mu_/R^, where mu„ — product of the 

gravitational constant and the mass of the moon (or sun), Rg*— 
modulus of the radius vector of the moon (or sun). 

The moin part Y^, containing the values X, Y, V^, Vy, V,.. 

or X, y, v^^, v^, v^; main part XS^, containing the values 

X°, yJj, or x°, y^, — directed cosines of the radius- 

vector of the moon (or sun); main part F^, containing the values 
X, Y, ”1, Vy, Vy, V„ or Xj y, a, v_, v„, v„, considering the 
other t.eiffls in the right members of enuntions (1.1) or (1.2). 

i. Results. In the main part F^ the new values of the 

right .lejuucrs of the equation/? of motion of the AES are addressed 
with the computation of the influence of Lruvltationol perturbations, 
produced by the moon (or sun). 

6. Algorithm: 

A, Vz ' (K-Z/r,)r,^/lr, W-zl) . 

h =/i + . i’y -Fy +As Vy , , 

lO -r\/rs = (d; -l/r, )' + ()'5 - Y/r,f , 


fg“~madulu£j of the radius vector of the moon (or sun), F^, F„ 

—components of other terms in the right members of equations (la) 
and (l.£), 

7. Text ^ 


mooUtlHt 
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00 I J«l*» 

' MOOf J>«VUI/R0 
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i 

i 


Oal./O/StOT*** 

?cj!>>*rCJ^3>*»005«lWCJ»*0-XO< 

XCTUOK 
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J) 
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1.6 Coraputatli-n of the Influence of the pres ure of light 
(F06-FL1GIIT) 

1. Designatian. The terms dolta^Vj^, deltaj^VY, delta^^V, , 

are determined, considering the influence of light pressure on tho 
right mevhers of equations (1.1) or (1.2). 

2. Gtruc.ure. Subprogram FLIGHT. 

General blocks: /CRS/j^. 

5. Conversion. CALL FLIGHT (KL, Y, RC, XS, KG, F). 

4. Initial data: KL— index, governing the computation of the 

Umbra of the earth (for KL=1 is produced the calculation of the 
influence of light pressure independently of the shadow, for KL= 

2 the shadow is considered); 

The main part Y^, containing the values X, Y, H, Vj^, Vy, Vr. 

or X, y, z, v^, v^, main part XS^, containing the directed 
cosines of the radius-vector of the sun: Xg, y^, sj or 

xg, yg, zgi considering the other terms in the right members of 
equations (1.1) or (1.2). 

5. Results. In main part F^ are addressed the new values of 

right members of equations of motion of the AES with calculation of 
the Influence of the pressure of light, 

6. Use of the domain COMMON, Before conversion to the 
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subprogram FLIGHT in the block COHMcN/BB/B it is necessary to address /16 
the value of the coefficient B (see par. 7). Starting from this, ^ — 
that in the kg, m, s system of units th^s coefficient has the 
dimensions m/s f it is necessary to mtke it agree with the system 
of units used in the coraputati^^ns, dividing by the scale factor 
from the block COHMCN/CEI.B/ELB (see O) , taole 2.3, N . 13). 

From the block COMMON/C one useds the constant (No. 10 table 
V 5 ) / • 

7, Algorithm: 


2*1 


''jT * • -X/rJ) r,Vl#’.-rl* + Fji , 

VV— a(y,‘’-y/r,)r,Vir.-ri*+/-y , 

Vg m-B (z z/r,)r,V |r«-n* +/*g , 
w. l.vrl/ra=^f-X/ref+^;"-y/ra)*+(zS-ZA^^^ , 


where — ^souare of the mid section, m — mass of the satellite, 

10“*^kg/m'^ — pressure of light in the region of the earth's orbit, 

kal for full optical reflection, k=1.44 for full diffused reflection. 

In the case of calculation of the influence of light pressure 
with calculation of the shadow of the earth the entry condition of 
the into the earth's shadow is verified; 

cts y < 0 M r rifi \ < R. t 

TX6 co&y -^ZZQ^/r , r , 

R — average radius of the earth. 
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8. Text 
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1,7 Right members of the system of equations of motion of the 
ABS. considering the normal gravitational field of the earth and 
standard five-la^jrardd atmosphere (FO6-FC1100) 

1, Designation, The subprogram FC 1100. is intended for use in 
the capacity of a subprogram for the right members for integratiun 
of a system of differential equations of motion of the AES 
models, characterized by the indexes 11100, 21100, 10100, 20100 
(aeo table 1.1), 

2, structure. Subprogram FC1100 

External subprograms: FNGRAV (FOl), RO(DOl), FATi:(F02), 

General blocks: /BHC/^ , /BK/r, /BSB/^. 

3, Conversion: CALL FC1100 (T, Y, F). 

A. Initial data: T— independent variable (time), Yg — main part, 

containing X, Y, Ya, Vj,, Vy, Vg, or x, y, z, v^^, v^, v^ 

5. Results: F^ — main part, containing X, Y, Z, Vvi V^, V«or 

X, y, z, Vy, v^,. 

6. use of the domain CCmMCK. In the block COMMCK /BSB/SB is is 
necessary to give the value of the ballistic coefficient, relative 

to the scale factor from the block COMMCR/CSSB/ESB (see (3) ), table 

2,3, No, 12). 

In the block COMMON/BK/KG, K6,KA, KS, KL there must be given 
the values of the indexes of tne model of forces (see par. 1.1), In 
the given case for the index KG the possible values are 1, 2; for KQ 
— 0, 1; the values of the remaining indexes are not used. 

In the block COKMOh/BHC/HC , as a result of the v/ork of the 
subprogram FC1100 is addressed the value IIG — elevation of the y^gg 
above the surface of the terrestrial ellipsoid. 


(J}u2a'Ud^ PAGK K 
op pt>>R QUA11TY 
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?. Text 


I 

f 
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l.T, Universal subOTogrsro for caluulatl.'n *»f right mrnhera of 
c<^uati’. nB a I moti'^n of *A1S with the use of various ;nodels of forces 
(FCT-fTRCE, Rt‘DKN3, LOGflOD) 

1, Designation. Subprogram FORCE is Intended for use as a sub- 
program of rirht-hand sides for integration of tho system of equations 
of motion of the AES for any models f forces, described in table 
1,1, The subprogram LUGIU D for the # iven model of forces determines 
tne value of the logical parameters, governing the computation of 
sidereal time, the positivins of the moon and sun. 

S, Structure, The package of subprograms, 

Ipputs: .. for the users: Fl’RCE, LOGMC-D, 

Internal inputs HCDEKS 
Utili55ed external subprograms: 


FNORAV(FOi ) , FATM {FOiy, Fa6rAV{F03) , FOR3S iFOA), 
FUOHT(FO^), AOJOACiBOb), RO(DOi), DENS{D02), 

SELESA (601), 37TCB05), OCLTLN(B\0)^ SUN (COVt) , ADEJJ , 
AMBAR, ORAV, TLOCAL{DOS)^ DE02, DE05, DEOAiCOi). 


General O^BiHe 0^OKiv<BA/i, /BlOO/^^ /BSB/[^ /BNAA/^ ^ /BRO/i, 

^T4, /vso4,/BDrm^ /BcmiR4^,/t«OEF/t^,/sYeAR4,. 
/CAED/^, /CAea/2,/aH24,/CBE/^y^Z4^A:AE/^,/CagLyi, /WSS/i, 
yCDSJ&'i , /C0MZ4 , /cauZR/^, , /'CAEl/i , 
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/COMy\, /CSDAY4, /0Ti4, /CP//J, /CEi/^, /BEM/^ ,/CeHO/(, 
/CHA/it, /BK4, /BrT/uj^ ^lECLyi,yMD/t. 

m mm * 


3. Ccm version CAL. Ft.I^CS (T,Y,F) 

•f. Initial data: T — Independent variable (time), y^^--main part 

of the functiuhs sought: X ,Y, ii, Vj^, V^, V-. or x, y, z, Vy, 


1), Kesulte: BV-‘*main part of functions derived from those sought; 

• ••• «v.^« •••. • • • 

X, y, A, Vy, V. or X, y, z v^^, Vy, v^— right me.iibeia of the 

system of difxereutiai equations (1,1) or (l.r), 
u. Use of the domain CCK.MCL, Before Conversion iu th* subprogram 
r’CHCI‘1 it iu necojisnry to ensure in the domain CiKr.rK .die valuos of 
Uie parameters utilised. In the block /BK/KG, KA, KC. K3, KL it 
is necessary to addresu the values of the Indexes of the chosen 
model of forces (oe- table 1.1), In the block CC-KintyBLOG/LAUK, 

LdUKil, ln5KL, L3!r by meons of Gonvereion to the subprogram LfGKvD 
(described below, in par. 3), it is necessary to address the valuos of 
the logical variableei, used for regulation of the computatit>n of 
sidereal time (for LHTvTRUE) , nosittons of the rai'on (for LGEL-TRUE), 
and sun (for LPUK:s'rRUE) . 

The remaining blocks CCKMCK are used only for the determined 
values of the Indexes of the .aodel of forces, for wnlch In each v.aae 
there is an indlcat hm. For KA giasater than or enual to 1 in the 
bloc . CcKM E/!'3B/GB It is necessary to address the value of the 
ballistic coefficient, relative to the scale factor from the block 
GimdV'USSB/^ (see (5), table F.3, No. 1?). 

For KL greater than C (computation of the pres ure of light) 
in the block it is necessary to replace the value of the coe- 

fficient B*(set‘ par* 1,6) relative tu the scale multiplier fr m block 
/GSLB/^ (see (i) ) , table 2, 5, Eo, 13) For KGsF or 3 it is necessary 

in the block COMKON/BNH/^to address the value of the number of 

harmonics considered in the gravitational fielvi of the earth, in the 
block GCLMC N/ BCOI’GR/, . ♦ —value- of the coefficients of the breakdown 

of the gravitational field of the earth (by conversion to the sub- 
program GGKGR(A02) (b) ). For KG::h in the block CCKEON/CA00/. it 


is necessary to transfer the values of the corresponding variables 
from the block /GA00A/2 (see (3)i table 2.1, Nos, 2, 5)» /20 


For KA greater than 1 or K3 greater than 0 or KL greater than 
0 inthe block CCKKON/BDT/^^ it is necessary to address the value of 


the date of tying the time in the mode RUB (relative to the Julian 
date). For KG=1, if KA is greater than \ or KS is greater tlr^n 0 


li 



or KL grente? than 0, and for KCsT, if KG ia grantor than , in 
tha block CfMMs i;/B3C/ GC, TG, KG it is nacasaary to addresjt the valua 
GO — 'BidaraaX tlraa kt midnight Greenwich of tha data DT, and TG and 
KG to sat enual to zero. 

For KAsT or 5 it is nectasory in tha blocks COMKOK/GKOSF/^q 

and /Gfh'AH/j to addrasa tiia values of tna coefficlanta of the model 

of the atmospiiere and set of nwarical corrocti ns for tha semi- 
annual aftect, which is acuorapllahed by conversion to the subprogram 
VKMACDOF); in the block one addres as the data and time 


tuc form of tiu? number of * h-hr. periods from the begimuat of 
^«'ar; In : i*e t lo^-k COMMDN/BRO/ilO? » FO, AP, D, 1 one &d ran ts 

of solar radio-radiation F,„,i wlvu 


In 

th.’ 

tuf value 


FI ;? 01 -^ho intonslt. 


computatiun of the 
leus thrn O.y), FC- 
( poo. ibis valuesj 
iii‘tic perturbation 
less than 0,b)» b— 


107 

retardation (for KA=h one nets F107 las 
average level of solar radio-radiatiun 

trihourly 


time of 
-of the 

100| IF:^, IbO), AP trihourly index of geomng- 
with cowputati. n of the time lag (for KA-7 set Ai 
purumetor, regulating ti»e calculaLion of Uie serai- 


tiian 


the 


at ion 


In 


- 0 ,- 

the 


annual eifeci (for 0 loo. than 0 the semiannual effect Is not con- 
sidered), I— parameter , regulating tiio calculation of the diurnal 
of root (CS) (for KC CS is cuneidered without tiie terra with ttie 
tiiO coefficient C , for I.-f CF is not considered, for I grrater 
0 CH in cunslderoa fully) For KA=4 (atmospher CISA-7F) in the 
block com K/BHchTlO?, FO, AKl , D, I it is suiMclent to ndures^ 
values FID7^ — ^Index F,^ n ittue of retardation 1,1? days), FO— values 
^107* ‘ivernged for fnbf' solar roiatii-ns, AKP--geomagnetlc Index 

considering ’>”“r,u67» GOO, 5^=3. ’'>33 add so on (time of retard- 

79 days), values of 13 and I without importance, 
block CtKMOK BHC/HC/HA0/HC, R1/B17CL/CS, GE/BFT/TM, 

GT, RG, AG, BG, RL, XG(3), XL(3), GG, GL In the process of operation 
of the subprogram Fi'KOF one addrea es the values HC— elevation of the 
AES abuve the surface of tha terrestrial ellipsoid, RC, Rl— radius- 
vector of AES and project I ms of the radius vector of the aES on the 
plane of the earth's equator (in the subprogram AGRAV) ; CE, 3E- 
Gosino and sine of the angle of Inclination of the plane of the 

of tii 0 ecliptic (In the subprogram 
GT — sidereal time: AS, BS— right 
the sun; R3, XG(3), Kb. XL(3)— are, 
and directed cosines (in tlie system 


of the angle 
earth equator to the pi ne 
TM=T--curi ent time (Moscow); 

■< consiun and declinatiun of 
ir. 'actively, for the modules 


UK) 


cu 

mu 


nates, determined by the index KC) 
GG, GL — mm-, / R5 and 

C(^MMON, enumerated In para, F, it is 
of the constants and sale factors in 


radius vector of 
In the remaining 

necessary to address 
conjunction with the 




the sun 
blocks 


of 

and 


the values 
tables 


'.1-F,h (!;) by means of conversiun to the subprogram GCKGT(AOl), 



7, The subsidiary subprogram RODENS, intended for Switching 
into the subprogram FORCE of the subprograms of computation ofthe 
density of the earth's atmosphere according to various models. For 
KGs:F or 3 one usee the model, realized in the subprogram DEKG(DOO), 
for KG=4 one uses the model CIRA-7F, subprogram ADEN (D03). 

Structure. Subprogram RODENS. 


16 


External subprogrami AQlQnC(B06), nCLTLK(BXO), DEUSCDCB), \PHn, 
*.MBAK, ORAV, TLrCAUi)C>5). 

General blocki; /BK/r^., /?>DT/^, /BR04 , AKOEF/q, /SYEAR^g , 

/BDVfAAR/j^, VCGM/j^, /CE^/^, /B»l/^, /CERO/^. 

U. nveralans CALL HcLELC (T, Y, IIC, RC, GT, AG, BG, 1). 
Initial data • r~“ti ®{ ■•"•X, Y , •<•, Vj^ , Vy* C R x , y , % , 

Vy, — modulus 01 radius vector and elevation of AIS 

above the surface of the terrestrial ellipsoid; JT- sidereal time; 
AS, BS«*rlt.ht ascension and decllnatlun of the sun. 

Result j p- -density of the atmosphere (in the uysteio of units 
used for the computation ). Concerning tiie use of blocks COKKtN 
see para, u. 


b. Gub| rogram LOGKOD, starting from the values of the indexes of the 
model of forces, given in the block G. MMCK/BK/KC, KG, KA, K >, KL, 
addresses In ti»e domain COKMrK/BLGG/LGUK, LGUKG, LGKL, LCT the values 
of the logical parameters, which may be used for the regulation of 
the caiculatltm of sidereal time (LGT) anu positi ns of the moon 
(L uEG) and sun (LfftlN). Below are introduced the Cimdltl?>ns, under 
whlcn each of these p.rameters takes on the value 


LGl'i.t WJn fur KG ,.ruAlfr then 1 or KL greater than 0 or KA 
greater than 1; 

I.aUNSrTRlLv foi KA greater th tn 1 or KL freAter th'-r C 
I.GHLk‘’'RUE for KB:rl ur K.Jr^; 

LAr-.RUK frsr K. ,reatir th. n 1 ot K-’, , renter than 0 or KL 
gre'^ter th”u 0, if KC=1 ana for KG grtster than T, if KC st. 


Uunveraiou: 


C a.L l OGW'D. 


8\ibroutini ro*CE<i*v»r> ' 
dimension V(S)tF<S>,VSO> 

LOGlCStASUNtLSUNS.LtELtCST 

COMXON/BrT/T»STtRS>ASfSS>lkf XS(S)fNL(S),SS|*W 

COMMOM/BlOSyLSUN«LSUNStLS EI.tiST 

COMMON/BK/KetSSfSAtKSvKl. I CALI, SUNl5T»X*ks»AS*SS*XSI 
* 


COMMON /COtt/ OlS»OOit 
COMMON /BHC/HC 

COMMON/bNM/NmO 
COMMON/cORL/OtL 
.common /BOT/OT 
IMT.EO.S) 

- OOTO 0 
T*X 

ir(iST) 

- . ST*STT(X) 

IF( .NOT.CSON) 

OOtO lo 


10 

10 


S 

i 


1 

ES 


8S*6llt/Rl/RS«0SRS 
SOTOftO,lO)»RC 
irusuNSi 

- call AO}BaC(ST‘»XS* tiXsl 
iri.NOt.LSiL) 

- OOTO 0 
CALL lELiNAlOT»X*XLiRL) 
6L>0 RL/RL/Rl 
eoTo<S(Oi*Kc 

call ACtSACf ST*XL(1,KL) 

continue 

call rN0RAVUCiK6*V*HC*R*r ) 
iriKO-ai iOliikiX 
CO >3 Jtli3 
V 0 («|)«V(^} 

60T0(2A,a5)»KC 

CALL AOIOaCCSTiVO* 1, VO) 

CALL rA6SAV(VtVe*r*Kf *N mO) 
ir(KA-Uii»E,U,|2 
CALL ROOSNS(X,V»HCiR»STiASfBSvP ) 
GOTO 3 

CALL R0(HC*R) 

CALL rATMUC»V*R*r ) 

ir.(KU10’*^03«13 

CALL rL!eHt(KL*V>R(Xf tRS*P) 

lrUS*l)l0A,Ei,32 

CALL rORSSf CSf V»XSiRS*r > 

eOTOl2l»iOA,3l),lCs 

CALL rORSSULf V»XL»RL*r) 

continue 

RETURN 

end 


I 



SyftMOUTlHt IIOOEMt(T»V«Ne*ll*ST««f*ftt*^| 

OIMlNSioN MUt«tO<A| lUXON|*> 

COMllON/ftK/KC«l(fti>*iKl*KL 

COMl^OM /ftRO/StOf 9 >'* 0 * I 

common /OOV/DT 

eONMON/BOVI«A/OTMi 


COMlION/ettM/OM 
COMMON /CI1/T5 


COMmON/cSOAV/ADAV 
COMmON/CC’/CI 
eOMNON/MM/eM»CtlC 
common /CCAO/EAO 
DO X 4m|*S 


I XfJl«Vfj)/A 

* 60T0(’«2>|KC 

1 CAU Ae|SAC(ST*X*2*X) 

2 SUI1>*A$ 

SUfCMBS 

satc9>»mc/cs»em 

ir(KA.C0*A) 

•. GOTO 7 

$ inou.sio 
9 0*Ot NC^t^ AOAV 

4 CALL 0CNA(SATlS>»X(l)*Xl2)»X|}|,Su*6CO|f»*tKOU»#Otl*OO) 

• A«Ao(«>*CAO 
GOTO 10 

7 CALLOClTl.N(X»SAt(2)»SATli)> 

AM40*D7*(T*T3 )/S0AV*i501*«G 

CALI. AoEMCAMJD*ty»SAT*CC0(TCM0«AC|0N*AMM7*»^} 

0«7/OM*iRO 


10 CONTlNUf 
ACTUON 
END 


subroutine LaCMoO 
C0MMON/.RLO6/4^5U|,.tSUNS,tStL»L»T 

LOCICAi. LSUn,i.$oNS#|.SEL»l5t 

CONMON/OK/BC,NB,KA*KS»Nt 

LOGlCAi, LSUN/^»t.r,»LStLR6 
COMMON/erT^T.STi 131 
T«»100 , 

lsunaska.ct.i 

tsUNs*ics.cT.|,o».ia,eT.« 

LSEL*KS.CO*1,ON,KS.EO«3 

Le*N6.6T.2 

LSUnblSUNA.OK.LSUNS 

LS*LSUN.OA*LtEL 

LKC«KC.CQ.t 

UST*LS.AN0.l.ir6.0*.LC. and. .HOTtLirO 

-Bimr*. 

INO 








Chuptfr Iiittgration of lysttai of dlfi#r#ntial touatlonn /2U 

(indax 0). 

fc.l* IritttrAfci^n of » lyntta of diffartntlal aquation* by th# .ath- 
od of Adaraa (QCl-AD/iKiJP) 

1. Dfol^notiaii* In tha aubprasraa tha aafchod of Adorns is raalixad 
with th# usa of d dlffsrancas, wltfi automatic choict of stap for 
ifttagretiun of * sy»st®» of dlffarantial aquotl-ns; 

Lrctj) 

with th# initisl conditions Y(t^)»YQ, whara Y, Y, Y^, F-ji-dimension- 

al vectors, t— 'indtpondent variable, Tha start is produced by the 
Runge-Kutta method cf the fourth order, 

?, atructura, Subprogram AUAHSP, 

On# uaae tha external subproitraas ; llilP and BKK— composed by the 
user, 

3. Conversion.* CAhh IDAKSF CKII, H, T, Y, im, BKK, KKY, 

MKY, El, B?, K, AF, FB, Fd, Tl, Yl, YF, Y3). 

/f. Initial data*. Kil-- variable, ensuring the possibility of 
Integrating with constant pact or with automatic cnolca of pact (In 
the first caae It is necas ary to take Kilal, in the second 
Kil=p); H—initial step of integration; T — first value of the 
independent vsrlahlo; 

Y— unidlmesialonal block of dimension N, contalninr the initial 
values of the deaii'od functiims; 

i Hi —none of the subprogram, computing the right mewbeia of equa- 
tions, set -up by the user; 

BKK--name of the subprogram, prepared by the user, intended for 
control of the end of integratiun. In this subprogram there occuro 
conversion from tho subprogram ADAMSl after each step of the inte- 
gration; 

NKY, MKY— number of the first and last functi n in .the block Y, /Eb 

respectively, according to which proceeds the control of th© 
accuracy of the Integr -tlon with automatic choice of step; 

El, Ec-blocks of dimension N, containing lower and up-er bounds 
of the per;nissible errors of the controlled f ractions; 

These values are used: BKNKYO, E2(NKY), E1(MY+1), E2(KKY=X)... 

El (MKY), EE (MKY); 

Ks.n dimension of the system of eauatlons 
AF, YE, Y3-* auxiliary Iglocks of dimension K: 

FS — ^two-dimensional block of dimension (N,S) 

Fq— twodimensional block of dimension (K,4)» 

Note 1, The designations of subprograms ERP and BKK must 
be written in EXTERNAL subprogram, used for subprogram ADAMSF, 

Bubprogravi i Ri must have the form Sb'BRt UTINK iRi(T,Y,F), 
where T — independent variable 

Y— block of current values of the desired functions 

F- -block of right members, length of blocks Y and F equal to N, 

Subprogram BKK; abBRuUTIKE BKK(T, Y, F8,IBKK, H), 
where — current value of Independent variable 

Y— cuiTsnt value of desired functions, block of dimension N, H— step 
of integration, 

F6— block of dimension (K,w), containing valuer, of ti.e rig. I member^ 


EO 


In ® cur §nt rotntfl nf intoii,r^»U n. 

I "KK— Inttrrjul wr.lcl;, within BKK, smut bt ticrlbtrf th» 

vala» £■ for fulflllfnint of t-.o condition*, Accofdinc to which i« 
dft#f«in#d th» momml of txit from tht lnt*gr«tlon (b«for* th.* 

•t*rt of Intifratlon th» aubprogram ^D/iMai ««crib#s to IBKK th# v/'X- 

ao I). 

In aubrrograa* iRi and BKK the vsl it* T, V, F , H ar» »itab- 
Xi«h*d by aubpro|sr«« At^AMiji In tu« proct«;« of inttgration, 

Rota 2, Jum* Uinta in tn# atarting portion thtr# arial* / 2 C 

a notU ir brtak down tn§ attp, not conittcttd with tnauring tht 
proacribtd accuracy of inttgratiun, Aacribini: to indtx IBKK tht 
vaiat 5 Itada to intorruptiun of tht proctas of atarting thnt 
was bift ,, dlsconntctlng th* influtnct of Kil on th# choice of atop, 
return to iht atarting point and to tht atart with a naived 
atep. 

Giving the index KBKK tht volat leada to a nr.v interruption 
of the procea- of integration, oatabliahmont of the initial atep, 
reeatablishwent of the influence of Kll on the choice of atep, re- 
turn to the Initial point and to the trnnafcr to the standard 
procoaa of integration. 

The subprogram Ai.1I3(HC4) may nerve as on illUHtration of the 
ni pllcftttun of tnie pusnibllity set par, 5,.|, 

b. Rosultiij 71 and the block values of the independent 
variable and the dsairod functions at the final point of the inte- 
gration. 

C, Hethod, /,Ve^ connider tho system (F.l), According to the 
known Values «d the Kth point and the ap;roxlmate values 

B^o^dced at the Kth and 7 preceading points, one determines 
by the extrapolation formula of Adame 

where h^-« step of Integration. 


Inteppolation formula of Adams 


si^ma 


aXluwa the obtainaU valuea of V be made moie prooiae. 

The coefi’icieiilB of the interpolation and extrapolation formulas 
of Adams have ti.e folloeiug valuea: 


ot, - 3,58995535, 
a, 4.9,5Zr:206Se, 
a, i 18,0845337, 
a, «>-22, 027753, 
a* * 17,mWoU, 
i-8,6I2I2797, 
a! « 2.44516369 


- 0,-304224537, 
p , « I,I5SI5S0b. 
p j a-l, 00691^4 , 

6 . « I,0X7964«Xt 
:M),7320S5S34, 

5 5 i 0,343080357. 

*^0.0938409392. 

Si ® a O.OTTRR7R942*. 
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For ©utoraatic choice of step of integration on© computes the /27 

# w— I- 

differences delta sYk+i — 

We give the allowable bounds of erx’or for integration: 

0 less thani or equal to epsilom ^ less than or equal to epsilon-, j. 

X I J f J 

If th®r 0 are fulfilled the conditions epsllon^^ ^ lea - than or 
equal to absolute value of delta Y^ less than or equal to 
epsilon^ j, where 1 ^ 2 — numbers of the first 

and last controlled functions, then the step of integration h^, does 
not change* 

If the absolute value of delta Yj is less than epsilonj^ ^ 

for all . .Np then there results a doubling of the step. 

For this Integration continues with step h^ until the accumulation 
of the necessary number of points, for which after an interval 
of time, the multlpTe 2h^ the known values of the function F, 

after which the step is doubled. 

If the absolute value of deltalTj^^^lj^ is greater than epsilong j 

even if for.- one value of j; less than or equal to j less than or 

equal to N 2 > then a breakdown of the step is produced. For this 

it is necessary to have the valu®*? of F in the seven (receding 
points with the interval h^/2. It is possible to obtain them by 

means of the Interpolation according to the formula of Lagrange 
for known values of F with step h,.. 

We introduce the variable xl =: (tj^ — t)/h., then the Interpo- 
latioji formula of Lagrange for the determination of P(t) Is 
written Jn the form: 

«t)»f*u7(fe-2l4‘+(7?l’-735l^+(624&*-l764S*+720t)/?O/.O- 
l -/i.6(^^-22^*+)90^’-820%^t849%*-2038%’+8^0^)/720 + 

. +#*.5(^-23t‘+270|’-925^''+2M.‘i^t2ii)2^‘'+(008g)/2«0 - 

-/■*.4(4-2'^^‘’+226^’-(056Ut25A?^®-29>2l|,*+l260^)/lJi« + 

+ -25^+24 7 ^’-12(9^'*+ 3112 ^ -3796 1 + 1680^)/(4^ - 

(4^-26 270 12, 20^^ 3929l^*-9272i g*+2>20 1)/2<0+ 
(|-27^®+295^’- 1665gk JlOA |.*-8028W+?0A0|,.)/720- 
28 £,*+322^^-19604^+6769^* 131S2^*+I3068|,-904Q)/50«). 


22 



To ^tain the valuas 
Runge—Ku^ta matsUod of tua 


of F in the first 7 points 
fourth order. 


one uses the 


/28 


£ + . , y*+i, ). 


l.itor/jtur *,*5 (n), p, 3. 


'If 

'W • 


7. Text. 


i002 

I 

lOOA 


1001 


100 

2 


'■p pt M g.- 


SUeoOUT iNE AOAMs^II^P'I »H» 2»V ,0|;K»MUV I tl » Ilf 

N*Ar»r»^i«xi*tfi»v<*v>i 

dimension 

Vl(N)*V2(N»*VStN)iA(tMtitt«0(*> 

«N«H 
|A»KPl 
H*MN 
lOiCKBl 

00 1001 J»i»e 
00 1001 S«1,N 
PU »J>«0 
ounce-kutta 
101 Z«X«0.T*H 
DO 2 J*1*N 

Xl’X 

.A(2)tA(ll 
A(3)«H 
0(4)*M 

AC9)«A(1) 

X2*Xl 
DOA Jci,N 
V2(j>sVltJ) 

-W)M J)»V1< J> 

00 19 lal»A . 

CALL PRP<**^»V2»Af> 

If* |•ll6»*#T 
60 toUO*T*T)«|A 
DO 12 JsltN 
rU.OlfAMJl 
If t lA-l)»#9rA 

CALL 6XKtXl»Vl»r»l»XK»WS^ , 

OOTOlO.t l000»1003tl002)»t£X4 

U*l/ 

H«H«,9 . , 

OOTO lOOA 
DO .0 X»l»T 

oa 0 . J*i' 


IP 

;r'-' 


6 

10 

12 


M)03 
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7 I) 

00 l9 j«l»N 

19 y2(J)>V9(J>«A(l)fArUI. 

tr(CXl*l)*Sl0N(l.iHm*l0iA4 
14 00T0a02*liil9)*|A 
14 tA«|A«i 
. 00 20 Jpl*N 

10 7l(J»l>BVlf 
.H»H«,9 
OOTO 100 

’ 19 00 17 JsNKV.mkV 

If uos(Vi(j)«ric^iiM«i2ijm7»it#ti 
17 CONTlNUI 

IA«l 

OOTO 101 
, 21 IA*|A-1 

. ooto 101 
102 tA*0 
C AOAMS 

1q 9 A(1>>-.30A224997«H 
At2)«2.449U309*H 
A(S)t-a.6l2i2797«H 
A(4)S17.379«544«H 
A(9>*-22*027793«H 
A(A)*16.0949Sa7«H 
A(7>«-9,92930a6a«H 
ACa)«3.9099S339*H 
“6fl>*;6il3673941*« 

0(2)«-.o93a4o9392«H 
a(3)«.343000397*H 
Of4)s-.732o393a4«H 
4(9>*l«ol29a4«l*H 
O(a)a-1«00491964*H 
0(7)81, i9«1990a*H 
0(0)8. 3 o 4224937«H 
' 104 CALI. OI(p(Xl,Vlf AF) 

00 23 jalrN 

V3(j)svKj) 

V2( jlavitj) 

CAilL*OKK(xl!!vl*r» tOKK(H) 

. 80To<109*1000il003»l002)tSOxK 

109 tr(|A-7)lia, 34*110 
34'H>H«2» 

00 39 K8l*3 

.L«0-X 

M 8 L-K 

00 39 j8l*N 
39 r(Jrt)sr(J»M) 

00 34 x«l»* 

• DO .34 jsliN 
34 P(J,K)8rl(J*Kt 

■00 90 JalfO . 

A(J)aA(j>«2, 

90 0(J)aO(j>*2, 

119 U*0 
lia Xl*Xl*H 

00 24 K»l»a 


00 't4 J«l*N 

•’iAiii oooUi«vt>An 

00 'lO 'K»l» ^ 

00 so J>l*N 

00 so Jal»N 

to J) 

ooTedooaoTitKst 

iOT 00 to JaNKV»MKV 

tf(ABS(vltJ>-Vt|J|U|tCj))tO*St,S0O 

to COMttNUf 

00 t# JaMitViMKV. ■ • 

II* Utg . 

Oofo 100 
so eONTINUI 
e HaMil 

in io-i)0&*oiiio0. 

. Si 00 09 K>1*0 
• k*t»K-l 
Do 90 J«1*N 
99 ri(J»K>>f(J*U 
e MaH 

100 00 90 Kal»T 

■ 00 90 jal»^N 

roe f u*K>ar(j«K*i) 

' tote 104 

C IfaH/a 

100 Xlaol.H 
laO.O 

Oo OT Kal»4 

■ MO)s-<ff(((Z«20<>*Z*922« | •Z-1*|0 « |•{•0740 • 1*1 

• *13 j92* )aZ4lO0Oa, )/S040,'^t*l • 

. OCTjaf I ( U { z-27, UZ«209i |•^•|OOo. )»t491o^» >*l 
-0020. Uz*5o4o. )/720««a ■ . ^ . 

. H(0>a-U U ( (Z-SO«)fz«27o« >•2-1410. laiaOttO* >*i 



•9274.)«z«2920.)/240taZ> ■ 

.||(9)«((f((l Z-29« )•Z•Z47t i•Z-121o• )f2«91l2« U| 
•3796. >«Z*l000. )/144«aZ 

; I ( ( (Z-24* }«Z4220* )• 2- 109*; >#242940 • >*2 

-2952. )«Z4l200u/l46«4< . . 

,.Rl3>«<f cn (z-23, >*Z#207«>*Z-020.>#2#2|444 >#2 

.. -24i 2. >*Z«loOO. )/240»*2 

. 0C2)a»(lt(((Z-22. > *Z ♦ 190 i > • 2 -*2 « . | #2* 1040 i >#Z 
-209*. )«z«040 . >/720.7Z 

. oa)«(C CU <2-21. M24l79j«2-739. |#z*1024, >#2 
•1764. )«Z«T20. >/5040f#Z 
00 30 jslfN 


n< j*K)so 

00 30lb1»6 
90 m< Jf K)*f 1< JtK>«l^<k>«f < J*l> 
97 2 bZ-1. 

DO 3O ac9,7 
L#K«K-0 


40 f(J»L)srl< J*K> 
Do 41 JsliN 

' V2(j)sy9<j) 

41 VKJ)*V3<J> 

DO 91 Jal.O 

4(Jt«A< il«.4 


DO 39 JclfN 
99 nj»t,>cr<J*K) 
* 00 40 Kc1»4 
ka2«K-l 


91 B<J)aO<j>*.9 
HaH«.9 
GOTO 119 
1000 RETURN 
END 


Interpolation according to Adams (G02-ADINT). 


n 


121 


1. l)©8ignation. The subprogram allows in txis process of integra- 
liun by the method of Adams of the system of differential equations 

YaF(t,y) , where Y, Y, F — -N-dimensional vectors, the computation of 

the values of the desired functions at the moments of time t^„, not 

nu’ 

multiples for the step of lntegrati* n. 

?. Structure, .jubpiurram ADIKT. 

3. Conversion: CALL ADIKT (A, F3, 7, YR, M, H, K). 

4, Initial data; 

A — val^e of the variable xla(tj^ — '^here t^^^ — given moment of 
tir.e, h^ — step of integration, tj^--current value of the variable of 
Integration, corresponding to the conditivui; t,^ — h^ less than t^^^ 
lesn than tj^; 

F8^ — block of dimension (N, 8), containing those produced from the 
desired functions in the eight last points; 


Y^-block of dimension N, containing the values of the desired func 
tlons at the point tj^; 

M — number of interpolated functions (M les. than or =N); 

— step of integration (h.); 

N — order of the system of^equations, 

3. Results: YR — block of dimension N, containing 


j~l , . • * ,M. 

6. Algorithm: 

(t 


v'/e Introduce the variable xis(tj^ — ^nu^'^'^f' 


nu' 


h^ and tj^ determined above). 


The value of the desired functions at the point t^^ is computed by 


the formula ^itJ 




where j=l,... , M less than or sN; f^^^ — values of the right members 
on the points (i=k-7, . . . ,k) ; yj(tj^)-valnes of the desired fun- 
ctions at the point t^| 
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V*-7 • + 700|®- 3528g®+ 9744^*- 1.4112^+ 8640^*)^ 

V*-0 — (2l|?-5284^+ 5330|®-275 s 4®+ 77658t^-II4I2^+ 70S6(4*J5 
V*-5 ■ C63^-I656|’’+I7388^-93240^4270I44£‘-4052I^+ 254016^ 
V*.4 — (I05^-8ae0£7'*3I640^-I77408|®-t«344S0|?-826560l? + 

> 4529200^), 

"VA-a “ (IOS^-aOOO|743458O^-2O4790g®.t65352O|?-IO6288O^ + 

4-705600^), 

'PA-2 — <63^-IW2fe’4a2680|®-I43I364®44950544*-886032|? + 

463S04Q£p), 

'Pa- I • (2ltf-648^’.^60^-55944^®+2I4368|?-449568fe®+423360&®), 
yf j, — (3^-96|'’’+I288|*-9408|®+406I4g‘-I0505^3+I568I6|? - 

-1209604) i 

[ 2 ] . 


subroutine A0INT<2ir,VfyR»M»H*N) 
dimension riNffiMVtNlrVRiNltRfS) 

Rm«UlUt3.»2-7 2.>«.2*7ep.I*Z-352S.)*E 

* ♦7744. . )*2*064o . >•? 

RI2)*((<n<-2l. •2*52 6. ) *2 -5 3 2 0 . j »2 4 27552 . ) •* 

- -77 658,)«2*-il4l26.>*2-7o56 0*U2 

l?3)«t ( I « t I 63. •2-1656. )• 2-* 17 386. U 2-73240. >#| 

* *270244, }*2. 605216. )*2*2540i 6« )»Z 
Rf4)i((((<(-io5. •2*2080. >*2-31640 «M2*17740B.>*2 

- -5S44S0.)^z*826560. )«2-5292oO , >*Z 

R(5)«l(((((io5.«Z-3ooo.|*2*3458n*)* 2-204 792. 1*2 

* *653320. ) *2. 1062860. > *2*70 5«00. ) *2 

R(6)«< ( ( ( ( (-63.*2*1672. >«Z*2266e.)*Z*14sl36. >«2 
. -495054. )*2*0B6a32. >*2*635040. )«Z . 

„ (7)*( M > < <21 .•2*640. ) *2*626 0. > *2 -35944. > *2 

* *214366. )»z«449566. ) *2*423360. )*Z 

fU6)«( 1(1 1 f-3. *2*96, )*Z- 1268. 1*7*9406, )• 2*40614. >*2 

* *105036, >«2. 156616, >*2*120960. 

DO 1 Jcl*M 

1 VR(j>sV(J> 

DO 2 Keli8 

R(K)«R(K>*H^2/12096o. 

00 2 J*l»M 

2 WR< J>=VR< J>-f( w,R>*R(K> 

RETURN 

END 


Chapter 3. SUBPROGRAMS, EI.'SURIPa IK THE PRCCSS • CB IKTSGRATIuK CF 
TliE SVSTSI-* OF CHDIKARY DIFFERENTIAL ECUATirFS THE ATTAINI'-Q OF THE 
PRB]3CRIBED VALUES OF THE FUNCTIONS FR' M THE SOLUTICK (INDEX H)« 

3,1, Arriving at the prescribed value of an arbitrary fun.*>,ion 
from the solution (HOI-REACH) in the process of integration, 

1, Designation, In the process of integration of the system of 
differential equatiuns (2,1) the subprogram REACH allows one to 
determine the value of the independent variable TR, corresponding to 
reaching the given value PRl of an arbitrary continuous function 
phi(T)sFREAClUT, Y) , Y — N-dimensional vector, 

2, Structure, Gubprogram REACH, Utilized external subprograms; 
ADIKT(602), subprograra“-functi.,n, FREhCH (composed by the user), 

3* Conversion: 

CALL RLhCH (T, Y, F6, H, E, FR, FH1,FHK, TR, FREaCH, K, YR) . 

A, Initial data; T — independent variable; 

Y^, — block of current values of the vector Y; 

iV 

Fa — two-dimensional block (K, 6) — -values of right members in 8 points 
with step H; 

H--value of the step of integration; 

FREACn--designation of the subprogram — function, computed 
phi(T)-- 

FR- -current value of the function FRSACH; 


*For the subprograms of this chapter the conversion must take place 
in the process of integratim of the system of differential equations 
i,e, from the subprograms BKK, to which is transferred the regulation 
from the subprogram ADAMSF (see par, 2,1) after each'%tep of inte- 
gration. 



FHK-^valae of FRSACH at the previous step of integration (at the 
point f — H) { 

FR1— given value of the function FH2ACKj 
K~ order of syr-tem of aquations ; 

v3— 'allowable relative error, miourlng the choice of value of the 
independent varlnble TR, oatlafy '« r.fj Ih© relationship 


FRi- f REACH (TR, YR) 
Af 


$ E 


wheiO delta F — Increment th the functim FiCH frcra the step of late 
gration. 

Note 1, Subprogram — function FRoACII, compnned by the user, 
must have the form; FnhMTi* F Fh.:ACn(T, Y) , where T — Independent 
variable; 

— block of current values of the vector Y. 
n 

Note r, Functlo FivFfiCH must be uef?Ci’lbed in the F\TFmFiL 
subprogram, unuig the subprogram RFACH. 

Note 3. Before. use^of subprogram RgACh it is necessary 
to find that moment of time in the pi'ocess of integration, for 
which is satisfiond one vif I'uc relationships 

F1?K loss than or -FRl less tt.an oi‘ -PR or FRK greater then or s 
FRl greater than or sFH, 

7 , Results: I'R- value of indeponuent varialle, ctu*responding to 

reaching the function FRFACH oi the value FRl; 

YRj^. — block of values of the desired functions Y at the point ?R. 

A, Algorithm. Let us consider, in the process of integration of 
tiie sy -tem of ecuations (?,1) a continuous function phi(T)=FREACH 
(T,Y0, where Y-'-K-uimenslonal vector, taking on the values FR at the 
current moment of time T, PRK — at the mo-ieni of time 

where H — step of integration. 

It is necessary to find th© moment of time, for v/hlch the fun- 
ction phi(T) takes on the value FRl, 

Under the condition, that thei'e is fulfilled one of the ine- 
Qualittes RK les than or equal to FRl less than or equal to FR or 
Ff?K greater than or =PR1 greater than or =FR, we determine F^rFR — 

FRK, xi^ = (FR— FRl ) /delta F. 

We set some epsilon greater than 0, If the absolute value of 
is les.: than or =epsilon, that moment of time TR=T — xi^H gives 

the solution of the problem proposed. In the opposite case, using 
the subprogram of Interpolation by Adams’ method (par. R.c), we de- 
termine the value YR at the moment of time TR and the value P-phi 
(1?R)=FREAGH(TR, YR). 



For th« new iteration w® deterrain® xl^-,+ (P-FRI )/dQlt« F, w© verify 

the condition absuluts value of xlple©-; than or aepsilon. If this 
condition is satisfied, thf?n TR=T — ^xi^H— xi^^H gives the solution to 
the proposed problem, oth* rwise we go on to a new iteration and so on, 

7. Text. 


»U»,0«TlHe "'‘“".I'ldiliiNi 

OimNUcN V(,)>rlN. 

opaf R*FRi 
V<«0 , 

Wl*0P^0Z 

call AO|NT<W.r»W'V«.N. »«*»*» 

FereK ACH* I* 

5p8p-f Rl 

60To’ 3 

RITORW 

ENO 
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3.2 L-atii-tainellun, in the processi of integration, of the minimum / 5fa 
and maximum of an arbitrary continuous function from the solution 
as functions of an independent variable 
(H02-EXT»!, iARAB) 


1, Designation. Dubpiogram SXTRM determines in the pro car,, of in- 
tegration of the system of equations (2.1) the minimum and maximum 
of an arbitrary continuous function phi (T)*FRHACH(T,Y) , 
where Y — N-diraonsional vector, and the value of the independent 
variable, corresponding to the minimal and maximal values of this 


function. 

?, Dtructure. Subprogram ESTRM, 

Internal inputs; iARAB, 

Utilised external subprograms; ADn:T(G02), FRENCH (composed by the 

U 0 0 • 

General blocks; /BMINA, /BFC/, , /BSM/.. 

3, Conversion; 

CALL EXTRM (T,Y,F,H,i1 ,1 2, FMlE,TMIi:,FF.AX,TKAX, B'R'r.CII, K, YR) 


Zf, Initial data: — Independent variable; 

Yj^ — block of cuiuent values uf vector Y; 

H—step of inlet ratlun; 

F— twodiraenulunal block (K, d), containing the values of ti*e right- 
hand sides at b points with step H; 

FI, i 2— permisb'iula relative errors, enaaring the choice of 
extremum, satisfying the ielationahips 


FIE-if(TE) 


AF 


<Pi 


vm 

OR 


TIE-TE 

H 


< P2, 


where B’IS, TIE- -value of the function and of the argument corres- 
ponding to the actual extremum; phi(TE, TS — approximate values of 
the extremum (if the functi and the coriesivndlng argument; deltaF — 
Increment In the function phi(T) for the step of integration; 

L — order of the system of eouati ons; 

FREACH — designation of the subpr^'gram- - f unc t i n, computed phi ( T) . 

'fills subprogram is composed by the user and has the following 

construction: 

rlTCTI. i: FRl’nCH (T,Y) 

whor<-‘ T--independent variable; 

Yjj— block of values >1’ the vector Y. 

Lote. roe name of the subprogram— fuuctl.m FMSiiCu must be 
doccrlbed in Ea foAKAL subprogram, used foi* subjircgram eJXTRM, 

5, Results; FMIN--last local minimum of function FhLiiCH for 

the interval of time from the beginning of integration t(.' tiie current 
moment; TMIR — moment of time, corresponding to the value FRSAGH= 

IKIR; FHAX — last local maximum of function FRSAGH in the some time 
interval, TKAX — moment of time, corresponding to the value FHSnCfe 
FMAX; YR|^ — block of values of vector Y at the point of the extremum. 

6. Use of the domain COKMC.ii; After the beginning of Integration 


in til® block CcMKi K/BKlK/Ht , IIB on® rauut substitute these vslussj 
HUO, K?sO. 

In tiite block CCKltrh'/BHC/j^ at each step of integration it is 

necaaaary to adureas th® current value of th® function FREACH, 

In the block C(HMCN/B3K/, the Hubprogram EXTRi' addresnea the 
vat a® of th® step of integral font 
7» Algorithm. 

Let HG, H?, HI be the valaes of the function FRKACH it three 
consecutive moments of time t^; t^st, tpst-h, t^st-h-h^ where 

h is the step of integratton, h^sh for integration with constant 

step (for integration with automatic citolce of step h, may be 


different from h)» 

If the Inequality (HE-HI ) (HC-H?) less than 0 is satisfied, ^ 
then the extremum (for HF-H1 greater tnan 0 — maximum, for HF-H1 
leso than 0— minimum) is found in the time Interval (t^, t^ ) . 

;Ve will term such a situation extremal. It remains to determine moi*e 
precisely the moment of time, corresponding to th© extremum, and 
th© value of the extremum. 

We proceed from the moraento of time to the moments tau^ 


= th.n tau,--0, tau-.l, taUj-l^h/h curro«,ond t. the 

values of the function HC, Hr, HI, 


ifvj. V CFii it ^ o.n jtuiL viJLav U»n L wiiib • ff& x viaw 


moments of time and the corresponding values of the function, 
starting from the fol i owing considerations. 


\MC-M2\<\H2-Hi\ THEN , S, • W2 , jr, « t, • > ; 
■Bt’.HC, ar,-t,- 0 ; Xj - I + yi,/ft . 

Ben \HO-H2\>\H2-Hl\nm t X,«t,»l8 

IF 

B,»KI, *j,= tj=l + A,//i ; Bi’=HC, *3=1, »0. 


Thus, we have 3 values of the funct on for three values of 
the arguments x^. We construct a parabola, pasoing through these 

points, and find the value of W, the argument of the extremum guaran- 
teud by this parabola. 

Turning to the subprogram of interpolation following Adams 
(par. r.c) with the obtained value of W, we determine the value YR 
of tn© vector Y at the moment in time TRst-Wh, then the value HR 
of the functlwm FREACH at this moment in time. 

If the absolute value of (HR- B^) /delta F is less than epsilon^^, 

or the absolute value of W-x^ is less than epsilon^, where 

epsilon^, epsilon^ — given values of the permissible errors, deltaP 

=HC-H2, then the problem is solved, otherwise we go on over to a new 
iteration, using the values B|=HR at the moments in 

time XjssXg, x^-Xi , x^sW and so forth. 
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m 


subhout iNi 

• rillACH*M*Vll) 

dimension VIN) 

OIMENIIqM X())*|0) 

COt'MON /ISM/S? 
t CCMI^ON /6MlN/HX#Ht 
common /orc/HC 

tf(H2U,2.i 

N2«hC-H2 

irlwl*lj2*»*»S.9 . 

♦ If tM;.i»*^<2.n*r<s,iM»o,i.,5# 

90 K*S|0Ntl* *Mi) 

ea»«H2 

|C9)«1«*S2/h 

Iff IHC-HX»«lO#»0»i 
0 0(2)«HC 
xc2)«e. 

5C5|*Ml 
ooTo xe 
f e(2)*Ni 

X(2|»X(9) 

6i9)«HC 

xOi*o 

10 call PARAMX>e*XM) 
f.«XN 

CAUL ACI^T<XM»r »ViVR»NtH*NI- 
tO«T-N«H 

HX«pREACH(ti,VO) 

tP((KX.biiM«K)29*2l»20 

20 If(AOS((KR*'B(ll )^l^'2 )-Pl|2l* 2t«2s 
29 .in AfiS( x^-xa ) ).f2)2l*21»2* 

2A SO 26 Jsl*2 
t« 4 - j 

XU)»XCl- 1 > 

26 e(L)sB(l'U 
Oll)*HR 
xinsx>/. 

CCTo lO 
25 .hR*B< 1 ) 

21 IF{K)32»S3‘33 

32 FMlN-HR 
TmIK*TR 
CCTC 9 

33 FMAXsHR 
TMAx«TR 

5 M1SH2 
2'h2=HC 
S2sh 
RETURN 
END 
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9. Tho subsidlery subproijram FARAB, lntind«d for th# dsttrmi- /kO 
nation of th« ixtremum of the parabola, pasting through threo pointss 
phl(X|), phi(xn), phl(x^) of some function phllx). 

Conversloni CALL IahaB (X, B, XK) , 

Initial data! Xj|— block ot values of the argument, block 
of valuta of th@ function. 

Hsaults: XH--valui ot th@ urgU'isnt, glvting tht ixtrtmum of 

tii® parabola, pas. Ing through the given 5 points, 

A‘*gor ithm, Let the valuua P, of aojr.t functi tn phi(x) be 
known for three valueo -f the argument Xi : -sphl(x.« ) . 

iJa> .#• mU 

U lo ntcctsury to determine .Oi'—tt.o valuo of the ergumtijt, givluf 
the extremum ot tne parabola, p&u.>lng thiough toe Indicated vaiueu 


m 




2. 

TaifiT* 




SUMOUTImC paSAB<K^I»KW) 

VsKtUtKlM 
DO Ijas*’’ 

TU-tUm J»-K* 

rbtuiin* 

INO 
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5,5* Exit of *u» Al§ »t ta« aHStnalnig node uf tht orbit (K05-K«. D*:) 


M 


1. binignutlon, ?uo sabprugrem permita one to determine In the 
pructaH of Integration of the s;^8t»«» of eouatlona U,X) or (I,?) 
the moment of time, Curronpondi ng to the -inr the aacetifHn^* 

node, and slno the coordlnatoii and components of ti.® velocity vaetoi 
of the AMS •it this moment in time, 

?0 structure, flubpioprsm NCliK, 

Utilized external subprograms i AUlUTCnr*') 


General dowalnj /B K/A,, /BltfA , 

9 n ...... , A 


5, C nvernloi 


CALI. r( .)?:(?, ¥,H,F8, TV, ?V, KB) 


if. Initial d'ltos T- time; 

— ^block of current values X, Y, .;, V.^, Vy, V. t 

ll--gtep of inteferutl !nf 

F 6 -»ttt 8 "aiin 0 itHioi.ni block of dimension (t,d), contalnlip the vaiues 
• ♦ # 

of aI Y, Vy, V. produced at 8 pointu; T-?U, ..., T-H, T; 

KB««loop numl r,''*; 

L, Heoultr.; If TV i» truator tnun 0, then TV, YV. —moment of 
time and block of v^^lues Xy, '*y» ^Xy corresponding 

to the A© crosiunp tne oHCsnuing node, 

L'oto, If the Guild itl n of exit of AIS ’it Ih© ascendlnf node 
lo fulfilled at the take-ofi stage by the raetlu.u of Runge-Kutts 
U,e, initial condlUunB «re given in tue vicinity of the nscen- 
di .g mnie), then the procedure of determining mor-’ precisely 
the mower. t cif exit at the node la disconnected, and thf 1 >op 
number (KB) is increar.ed by 1 * 

f>. Hue of ti.o dn.ualn GOmu, 

’Qfore the beginning of Integroti n in tie block C *" K/B’,k/''K, 

LKL it Is necassnry to address the value Ht ) to the place -K and 
IKK; in the block M-u IV’3IN/IN0D on^s sets 
V. Algorithm, 

In the pj’ueesB of integration of ©nuat’ ns of motion of the Att 
one deteiiunes ‘he moment of time t, for whieii Is fulfilled the 
conditions A(t-h) ’.(t) les than U, ".(t) gronter than 0 , 

.Itartint fro.m the value of the parameter xi, :'.(t)/delt;« :.5 
whoi e delta W/.(t) — X(t-h), using t..* .Subprogram ' :»f inte! i iL ; ui 
acuurdlng to Adams (par. B.B), we uelermrne the values Xy, Yy, :iy 


Ay 


Vtr , V.. at the mument of ti.ue =t-xih, wiiere h — stop of 

X y «Jy li\X 


Integration, 

If the absolute value of 3y/delta 3 is less than epsilon as 


given, then ty, Xy, Yy, Xy, Vj^ , vy , 


Vy , V., give the solution .0 

'■y *v *'v 

the rrnnosed problem. If th« conditl n is not fulfilled, then a 


new iteration is carried out, starting from the value xl^,“xi 


t 


rt- 





.'iy/Uolta and so forth 


Note, Th© chossn value epsllonslO"'^ (variable BN, value 
of which Is given in the operator DATA) ensures exit at the ascending 
node with an error, nut excee'iing 0.01 m. for the ASS with apogee 
on the order jf 100 km, 

o. Text, 


SUAKOUTINC N00ifX>V>H«r I tv »VV 
0|l!)ENttON VUI*VV«AI*r(«*AI 

eoMMON/t|K/|K»IKN 
eOMjAON/etN^tNOO 
rvsk|o0»- 

KsiiaHdtfN) 

MsvmetK' 
u irMUji#*'- 

1 irMsstrci*i»*rc»#iwr< 

.10 OOTO. 

It iNODst-* 


loi 


N0>NB41i 

•ofoi: 

T OI*A»>ntE* 

M«lfCn/0< 

f etkt< AOiNTtN’fr*v*vv*t«H*ti 
>I1V¥VCB)/0I. 

J. . , : 

irP ABSiwi>^fKiiiiii*o 

U 

» zR"vor« 

KifUIUI^ 

END' 


n,‘f. uetci';..M»ution In the process of integration of the minimal and A’ 
maximal elevation of the AES above the surface of the terrestrial 
ellipsoid (HG4— AP3IG). 

1. Purpose. The subprogram APdld allov/s the determination of 
the minimal and maximal value of the ele vation of the AES above 
the surface of the terrestrial ellipsoid in the proces.j of integra- 
tion of the eauations of motion of the kl^ (1.1) or .(l.P), 

2. Structure. Subprogram APSIS. 

External subprograms employed: HSIGHT(BlO) 

ADINT(G02), PARAB(H02), 


■f- 
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Common blocks; /XKA, /B^lS/y /BKC/^, /BSB/^, 

3. Conversion: CaLL AiClC (T, Y, H, F6, HMIN, UMAX, 1BKK).* 

k* Initial data; T — time} 

Y — block of current values of X, Y , X, V^, Vy, V- ; 

H — stop of integration; 

FS-- two-dimensional block (u, 8), containing the values of 

X, Y, X, Vjj, Vy, A derived at the momonto of time; T-7H, T-uH, 

BKK — variable, controlling the exit from the integration by the 
subirogram IDAMdl (par. F.l) 

5. Rf suits: 

HMIN, HMAX-->minimal and maximal values of the elevation of the 
AES above the surface of the earth's ellipsoid in the time Interval 
fr-v- the beginning nf Integration (or fmm the start of the loop) 
to the cur ont moment of lime, 

Remark 1. In order that HMIK, HKAX be determined as minimal and 
maximal values of the elevations not on the whole interval of 
integration, but ou each loop sep a*aiely, it is necessary to address 
the values HC— of the current altitude of AES ' -at the moments of 
exit of the ASS at the orbital node for HMIh, -ItXAX, 

t.. If the extremal situati.un arises in the starting 
portion, then the varinule IBKK receives the value 5. 

in this case tuo program of integration ADhMIP automatically 
reuuces the step in half, switches out the influence of KPI on the 
choice of step and begins the start from the initial point. After 
HMIN or HKAX is found, IBKK receives itie value if, the initial step 
of integration and the influcence of KFI on choice of step are 
automatically renewed, ana a new st'^rt will begin from the initial 
point and the integration enters the standard regime. 

6. Use of the domain CCKHOK. Before the beginning of inte- 

gration in the blocks CCr’IOK/BSK/.'lK, XKKyBAl 3/1 HK, HI, H2 it is 
necessary to give the value X(t„) to the variaules ZK, XKN, to set 
IHMsl, Hlt=0, H2=0. ° 

In the block CrPKOK/BHC/HC at each step of integration it is 
necessary to address the value IIC, current elevations of AES above 
the surface of the earth. 

In thf? block C0KM(-N/B32/, , the subprogram A13IS adureases the 

value of the step of integration. In the block COM'‘-XIv:/BTAPS/ 
TMIN/TMAX are addressed the moments of time, corresi'ondlng to 
HKIK, HKAX. 

7. The Algorithm is completely analogous to the algorithm, used 
in the subprogram EXTRM(H02), only in place of the function FKFACH 
one uses the altitude of aes above the surface of the earth's ellip- 
soid. 

For the permis ible erx'ors e^re chosen the values epsilon.j = 
0.0001, epsllon^^O.Ol (variables PI, P2, values are given to tnese 
by the operator DATA), 
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subkoutinc •^*^t***M*****<***^>. 

OIMINStON 
OINfHBlOH V|6) 

DiTi 

eoMHON/BSt/|B . - 

COHMON VBAyt/tMM*N|,M 
COMMON/BZK/tlKltKN 
eO'MMON/BMC/NC 
eBMMON/eTB^S/TM|fi*TMAX 

u NMtNsHC 
NMAXaHC 
• OTO t 

1 tr(H|l9«9*’ 

3 

W2«MC*H2 

trCN&*»(f »4*9*9 

4 irirc >2o»3i»3o 

.31 60TO|91*9! r3»*|MW 

91 tMMhS 
lftKK«3 
IB Hl«0 
H2«0 
2KS{XN 
BETUfN 


3«; XBSlCN(l.*»<lt 
6fl>«M2 

X(3)sI,*S2/H 

irc ( hc»hi)«k)9*b«b 

• B<2)«HC 
XC2)sO, 

9(3>«M1 
60T0 19 
9 B(2)-«H1 
X(2>sX(3) 

BC3>«MC 

X*3>»0 


10 call 9AfiAB(X>e.XM) 

call AriHT<X»’»f »V*9 R»**Hi6) 

CALL HE|CHT(VR ,hR > 
ir <(UP.R< 1 I )«K)29,2l,20 

20 ir(ABSC(HP-B(in<'»'2)-Rl>2l*2X»23 

23 ir<Al!S(XM-X(ll).P2)21*21i24 

24 DO 26 Jsl*2 
L*4-J 

XfL)«X(L-i) 

26 B(L>sB(L-n 
B(1 )rMP 

X(1 )sXM 
GOTO 10 
29 HR>6(i) 

21 ir CK)27,28»28 

27 irtHP-MMlNI32»29»29 

32 HMiMsHP 
TM1 nsT-XM*M 
GOTO 29 

28 If fHr-HMAX»29»29»33 

33 HMAxsHP 

TMAxsT-XM«H 

29 GOTo(9>93*9)tlHM 
93 IMMs7 

IBKKsA 
GOTO 18> 

9 H1SM2 
2 H29HC 
S2«H 
RETUPN 
ENO 


CHAPfKfi 4. 


JUBPRCGHAK3 FOR 
FUKCTI' KG FRv K ' 
C0X 1). 


OOKPuTiiTiOK vF 
” ; i ARAMSTiRG 


THG VALUKG OF DIFi* FRONT M 
OF MOTION OF TUF ABS (in- 


4.1. Right ascension of the AES , local time and zenith distance of ^ * 
the sun (101— RIGKTA, TLCC, F.DSUNL, nDGtINC). 

1, Purpose, The subprogram TLCC for known geographical longi- 
tude (lambda) of ABS and sidereal time (GT) determines the right 
ascension of the ABS . 

The subprogram TLOG for known geographical longitude of the 
point (lambda), sidereal time (ST), end right ascension of the 
sun (alphag) determines local tine. 

The subprogram "DGUKL for known geographical latitude (phi) and lon- 
gitude (lambda) of a point and for the directed cosine of the radius 
vector of the Sun (xg, y^, Zg ) in the Greenwich system of coordin- 


ates detei'raines the zenith 
The subprogram 2D3UKC 
for known directed cosines 
earth's ellipsoid (x^., y^, 

ius-vector of the sun (xg, 


distance of the sun. 

determines the zenith distance of the sun 

of tne normal to the surface of the 

Zk ) and for directed cosines of the rad- 

y^, z°)in the Greenwich system of coor- 


dinates. 

2, Structure, The package of independent subprogram- functions: 
HI GUTA, TLOC, GDGUKL, GDSUKG. 

Common domain: /CVl/y /CDFGR/, , /CHRAD/^, 


Remark 1, Gubprogram — function GDSUKC refers to tlie subpiogram- 
f u nc t L 0 n GDGUKL . 

3. Conversion: A=RIuIITA (ST, ALK), 

TLcTLOC (ST, ALN, AG), 

SS=SDSUKL (ALT, ALK, X3G) , 

SSsSDSUKC (XKG, XSG) . 


4. Initial data; ST — ^sidereal time in radians; ALT, ALK — 

geographical latitude and longitude of the point in radians; 
AS--right ascension of the sun in radians. 

Block XKG, — directed cosines of the normal to the surface of 
the earth's ellipsoid at the point underneath the satellite; 

Block XSG^--directed cosines of the radius vector of the sun. 

5. Results: A— oocensiji-, of the ABS in degrees; 

TL--lucal time in hours; 

MG— zenith distance of the san in degrees, 

f?. Use of the domain COMMON, The constants from these blocks 
are used: COMMON/GPI/, , /CDEGR/^ , /CHRAD/^ (see 930, No. 1, 2* 

3 table 2.3). “ * 

7* illgox'ithm. Right ascension of the ABS,— lambda , local time 
— tg and the zenith distance of the sun — hGg are determined ac- 
cording to the following formulas; 
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* 


» X +(S7' , 

« A* + JT 

C08 ar J + ^ ar ^ . 


3. Texts. 


M 


function R|CHT«CST(*i.N) 

COMMON /COECR/OEOft 

common /CPI/mI *PID2*P12 

T«St*ALN 

R«T/P|2 

T*T*R*Pi2 

R16hTA*T*0EgR 

return 

END 


! FUNCTION TLOCIST* ALN*RS) 
C0mm0N/cPI^P1»P|02*PX2 
COMMON /CHRAO/HRAO 
TsSt*ALN"AS«P1 
A IF(T>A,5»9 

4 TsT»pi2 
POTO A 

5 ksT/PjZ 
T*1-K«Pf2 

tloc»t»h**ao 

RETURN 

END 


PUNeTfON 20 SUNC(XE»xSC) 
OIMENilON XE(3I*XS6(3) 
COMMON/CPl/PI»Rl02iPia 
common /COEOR/OERR 
c*o 

DO 3 j*l*9 
J C«X|I J)«XS8( J>»c' 
S«S0RTU-CRCI 
. lr(C>7«R*7 
R ZS*pl02 
60To. 11 

r ZS*ATANf S/C), 

irczs)io(ii*ii . 

10 ZS>ZS«P1 

11 ZDSUNC*ZE*I>C0R 

return 

eno 


FUNCTION ZDSUNt(AtT»ALN»XSO) 

dimension XSOO) i^XRIS) 

C8«C0SlAtT> 

XR(1)sCoS(ALN>«CS 
XR(2)«SiN(ALN>*0D 
xr(3»»sin< Alt* 
Z0SUNLxZ0SUncCXR*X$6) 

RETURN 

eno 



Geomagnetic coordinates of ti^e AES. 
TGEOM). ’ 


geomagnetic local time 
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1. turpoat. For iu» point, given by directed coalnee (Xj^, yj^, 
z,t) in the Greenwich oystara of coordinates foi the normalto the 

earth's ellipsoid, one determines the geomagnetic latitude and lon- 
gitude (subprogram GKLLO and geomagnetic time (subi-rogram TGECM). 


r, structure. Subprogram GMLL and subprogram- funct I on TGECM, 

General domain; /CII/3 , /CDEGR/^ , /EGU./^, /BGS/.,. 

5, Conversion; C.LI. GHLL (XRG, GLT, GIF), 

TKsTGECK (XI:G, X3G). 

q. Initial data; Block XXG,, — ^dlrected cosin^c of the normal 

to the Surface of tho earth's 0liips,?L<i In the Greenwich syste.;. 
of c;ordi nates; the block XBG? — directed cosines of tue radius vec 

tor oi tne sun in the CJreenwicu system of coorainntes. 

y. Results; GLT, GLN- -geomagnetic latitude and longitude in 
degre- s; 

TM — local geomagnetic time in hours. 

o. Uso of *viio domain COMKiN. From tho blocks Ct’Kt'.OX/Ci'I/^, 

/CdEuH/^ , the constanin a e used (see (>), Kus, 1, r, table £.3) 


In ti.e block CCkkuR/BGLL/GLT, GLK the subprogram TGEOK 
gives the value of the geomagnetic latitude and longitude of the 
point, in the blocs CuKKuN/HGG/GLT, GLR — the values of the geomag- 
netic latitude and longitude of the sun (in degree's). 

7. Algorithm. The geomagnetic latitude phi^^ and longitude 

larabdCjjj are computed by the formulas; 


. 0 


te 


.0 

;« 


.0 

:8 


•xiiXg +yjy| ■*s:S*S 


yzxi -X7,yz 


0,7I44707B. — O.I86I2600X,i;^ >>0,979924705 - 


m 


directed cosines (in the Greenwich system of coordinates) of the 
radius vector of the geomagnetic pole (phiQ=73.3°» lambdaQ~-69°) . 

Geomagnetic time TM=lambdajj^^ lambda^j^g, where larabda^g^ — geomag- 

netic longitude of the sun, calculated by the same formulas, as 
lambda^jj, if instead of y^, one substitutes the values 

Xg, y^, Zg, directed cosines of the radius vector of the sun inthe 

Greenwich system of coordinates. 

8, Texts. 
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subroutine 6Mkk(XitakT»oW> 

DIMENSION XE(S) 

COMMON/eOEtR/OlSR 

er«soRTU-»r*sri 

8bf«S|SNiR(Ot*SF>«OSSR 

sofo 5 - ■ ■' "r' 

.SkT«4tRNiSP/er MOSER 
MsSr^.fTftSiTOS.REtS) ^ , 

er*XE<2>>tETt44Te7S«.ld«il2E0ttl»M®t^) . 
lriMlS«4*S 
abN«SlGNCRl02tCr> 


OOtO • 

S BLNsATANUF/W) 
IF(|V)7»B(E 
7 ' 6 LNb 0 LN*F 1 
GOTO E 

A IF(EtN)9*A«E 
• SlNiSLN*Fll 
A OLN«GI.N«OESR 
RETURN 
' END. 


J 


1 

2 


f V 11 ^ 


• V 


COMMON /BeS/SLTfSLN 

COmmON/bGLI/cLT > OLN 

dimension KE(S)iXS6CS) 

call OMLL<XEiCLT»CLN) 

CALL 6MLL<XS6.:$LT#SlN) 

W»(6LNrSLN)/15i+l2, 

H'(N)lf 2 f 2 

W«2A,4W 


TCEOMsw 

return 

END 


.ffflSA|SS«XSIS} 


If. 3* Auroral longitude, auroral time 

(103— AULONG, TAUR) 

, 1. iurpose. For the point, given by directed cosines of the 

normal to the surface of the earth's ellipsoid (xg, y^,, ) 

In the Greenwich system of coordinates, the auroral longitude 
(subprogram AULONG) and auroral time (subprogeam TAUR) are determined. 

?. Gtructuro, lJubprograra— functions AULONG, TAUR. 

Common blocks; /CPX/y /CHRAD/^ , /BaU/I. 

3. Conversion; AUsAULONG (XNG), 

TAsTAUR (XNG, XSG) . 

Remark 1. Gubprogram^ — ^TAUR refers to the subprogram AULONG. 

4. Initial data; The block XNG 7 — directed cosines of the nor- 
mal to the surface of the earth's ellipsoid in the Greenwich system of 
coordinates; block XGG 7 --directed cosines of the radius vector of 

the sun intbe Greenwich-^ system of coordinates, 

S*. Result--: AU-'-auroral longitude inhours; TA--auroral time In 

'^L Use of the domain GOKMCN. From the blocks COMltON/CII/,, , 

/CHRAD/^ ono uses the constants (see ( 3 ), Nos, 1, 3 table 2.3). 

In the block Ct;MM( K/3AU/AU the subprogram TAUR addresses the 
value of the auroral iungitade, 

?. Algorithm, The auroral longitude ^dcording to /31 

(b) is computed by the formulas; Ka - anstg (A/B)-56*, 61 , 


. = 0,930414 .Vg - 0,331032 - 0,167062 , I 

= -0,353629 Xa - 0,923585 y» - 0,147963 . , ; 

^ a x% ^ 0,062724, 

i » 1 0,022319, _ , 

a » -0-,0I0589. 

Auroral time g 

computed by the same formulas as instead of 

x^, y^, z;^one substitutes Xg, y|, Zg, directed cosines of the 

radius vector of the sun in the Greenwich system of coordinates. 

8 . Texts. 


page t!f 

/'j'- 
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•rUNcTtON AUU0W6<*«* 

C0M>4eK/cHRAD/Mt*0 

* OIMCNSION XfO) 
xitxe* n*.8*iTt4 
)(2*)(E(2)'”t882Si9 

X5*XEt>l*.0l8SJ5, 

A»,9SC4i4*Xl*»S8lCi2»X2- 

B«-.S5J*2»*XI-.8»3»«5*K2 

2 M«S|6Nlpt02»Al 


i9T0*2«XS 
,14724 j*kJ. . 


aOTO 9 

I MsAt4N(A/2) 

IMB>4,9»5 

't, 

60TC 3 

9 IFCh)6;3i3^ " 

*6 W ^ H ^ 1 2 

3 AgtO«^O«{W-t98e030e9)*HRA0 
RltU^N 
2ND 


rUNCTlON TAU2(XC*-2ta) 

OIMtNStoN XE(3I*XSB(3> 

COMAlON/BAUi^AU 

AU>aUI.ON*(X|) 

AUS>AULoN8(XSS> 

M«AU«AUS^IS* 

Xr(MU*2*2 
•l W>24»^W 
2 TAUR«W 
R8TU2N 
END 


if.i't, Goomagnetic parameters B, L (104- BL, INVAR, LINHIG, STAR, A 
CAK-iEL, NEA'MAG, IKTEG, ASIN) . 

1. Purpose, Subprograms intended for the calculation of B, L, 
the parameters of the magnetic field of theearth, 

2. Structure, Package of subprograms, Inpuj: for user: BE, 

Internal entries: INVAR, LIKES, STAR, CARKEL, NEV/MAG, INTSG, ASIN, 

Common blocks: /P1/^ , /KAG/^204» /COP/y 

3. Conversion: CALL BL (XC, YC, AG, NG1 , TE1 , B3, FL), 

4« Initial data: XC, YC, ZC — Greenwich coordinates of the point 

in meters, KG1— number of harmonics considered (NG1 lesn than or = 

9), TE1 — period, for which aru computed the coefficients of the 
expansion for the magnetic field of the earth, minus I960 Uor 
example, if the period 1963 is chosen, then TEI:s3). 

3 , Results B3--block, containing four actual magnitudes: 

^phi» ®theta» ®r» coordinates of the vector B In the geographic-il 




nystem oi cuiirdinattib and tha mudulUB of th® vector B in dat>r® S| 
i?'L--V'fiae uf L In rodii of tt.« snrth, 

6* Us® of domain COHMoF* 

In tiiQblock CCWlObyCDB/j ar« adursa ®d 11^,, phi^, lambda^ -* 

geographic coordinates of the magnetic — conjugate point (Hq in km, 
phiQ, lambdaQ— in dai.rae8), 

7, There la a more detailed descrlpti'm in (6). 


t+.&. Ittvaiinnt geomagnetic latitude, 

( 105 — alKLAT, 01KL‘T) 

1. iarpoae. The subprogram AINLAT determines the invariant 
geomagnetic latitude for known values B and L, of the geomagnetic 
parameters. The subprogram OINLAT determines the invar Ian f latitude 
for the basis oi tiie line of force. 


P, structure. Subprograms — functions AINLAT, ciNliAT, 
Common block /CDEGK/^ . 


3. Conversions: 

4, Initial uu<-a; 
field of the earth. 

Results ; ALI- 


ALIr/ilKLAT (B, FL) , 
ALOsCIKhAT (FL), 

B, FL — B, L- parameters 


!l the magnetic 


liLO* 


^ , 


'invariant 
o. Use of 
block /CDHGR/'^ 


IcVtituue OI 


invariant geomagnetic latitude 

U \f% r* ^ #x rx 4* -}' 

liit* domain COMMON: One uses the Cimstant 

(see (y) No, T table 1,3) 


indegroes, 

wXVv X.11 

from the 


7, Algorithm* Tiie invariant geomagnetic longitude Lambda, 
according to (C), in computed according to the formula: 


/V ” arctg( ^i/Bg "I ) > 

rae Bg 


. 1 , 259921 , 
>^, 1964259 , 
b^ ^, 04686632 , 
bg ^, 01314096 , 


4 , . 4 , 00308824 , 
b^ m 0 , 00082777 , 
6 , .- 0 , 00105877 , 
bj . 0 , 00183142 , 


(Q.an653/B)M3 

L 


HaBapnaaTHaH nuq>oTa y ocaoBaaBa caaoBott aiBBB, 
Aq * arctg (/T=l) , 
rae B,L - reoMaraaTiaie napaiieTpH. 


S* T®xfcs» 


ruNcTioii oiNiATini 
COMMON /eDECR/DltR 
01NIAT«6 

lrCTN)t*l*l 

t O|NlAT«ATAN<SOllT(TNn«OE 0 ll 
2 "ITuEh 
end 
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PUWnTlOw AIWt.At(!R»Fl,> 

common /CPEftt/OSO*t 

A«< .3Xl093/0>««'l.3S333S3J3>/rl. 

B5s*.U<(fl CA*.fl01«Sl*2-.OOia»«TT»#A*,000**TT7)*A?.B038B»*‘» 

« •A-,8l3l*C*a>**-.0**Na*»*>**-.‘’®*29*****>^***’**** 

AINtATRB. ^ .• 

?M# j ,/ea«i» 

It I? ( T N ) 2 . 2 » l 

A'lNl.AT*ATANC*ONT«TNn»Of «N 
2f TURN . 

END 


*1,6. Geographical coordinates of T--ioint — points of inter- /bk 
section with the plane of the terminator of the line, connecting 
the Aiss ’vith the sun ( 106--FTERM) 

1, i’urpoee. From the known Greenwich coordinotes of theAES 
(x, y, x) and the directed cosines of the rSdius vector of the sun 

(Xg, y^, ) one determines the geographical coordinates of T — the 

poioit, 

B, Structure; Subprogram FTSRK, One uses the external sub- 
program GEOGRC (BIO). 

3* Conversion; CALL FTERH (XG, XSG, I!T, TLT, ?LN). 

If, Initial data; Block XG^ — Gre nwich coordinates of AES ; 

block XSG— directed cosines of radius vector of the sun in Greenwich 
system of coordinates. 

5. Results: HT, TLT, TLN- -elevation above the surface of the 
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aarthj latltaae nnd lonfituu® of t In radtann). 

f), Alporithft! Xr^t z..— Oreinwlch coordln.-ttes of the T- 

point, points inttrisction ,vltu tne plans of ihs tarmlnotor of 
tha lint, sonnasting tha 15.0 aun (or co:.tlnuatlon of tula 

Una), art computed with the fonnuiaa; 

f y-atVl, Zr^K-at*l,' 
a^mxxl *yyl +**V i 

The gtografihlcni coorainwtoo or T--point » • .• •• • : i. ih© 

sorili.rn par, 3,10 (I?), 

7, Text, 


5U5R0UTINE PT»- • ( XC » X*6 • JtN) 
DlMCNSlo*^ Xq(3><«XS6(S)»XT^ — 

DO 3 Jeii5 
3 'e*C4X0( j j i 

Do A j”! * S 

4 XT{ J)sXCU1»C*XSSU) 

OAWL OEO&RC(XTfHTUl.T *Tt*<*XT| 

RETURN 

end 


h»7, D«tar‘;ilnnt I . in f the ^Irer-nwlcii co"f *li ivitar; of tiis point, 
given gsographlcnl Intitude, longUudo ?md el«'Vf«ti 'n n^^ove the ear- 
face of ih enrti.’s elilpouid, and matrices of ti.e trnnsformnthm 
fri-o tho nreenwich to the; point sy. tern of coo, din ten, connected 
with Inis point (I 07-*h;lhI'), 

1. i'urpoo.e. The Gre nwlch uaorulnaLeo yiG» ?’'j.g)are 

datonHined for the point 1, pi. ascribed geof.rn|i;Lcal latitude, (phip 
longitule (lambda^) and elevutivn (h|j) above vue surface of the 
earth's eliipeold, 

file poiiii topocentric system ta" coordinates, connected with the pulut 
i , is determined in tne follow!, g manner. 

The origin of the cooruinoies coinciues with t«e p.oint 1, tno 
axis i X. is directed to the ni>rtn pole oi t!io oartii by the tangent 
to the meridian of the point 1, U.o axis I’yj. "-hy the external normai 

to the surface of tiio earth's ellipsoid, ;.ho .-.xls i Zp completes itte 

system up to tue rignt aide, 

Tne subprogram iOIMT determines the matrix of the transformation 
from tue Oreenwich system of coordinates to tifO point system. 

2, Structure, Subprogram 
Com on domain t /CAK/-. 

f 


5. Gonversiunj (lijr, it-T, iLN’, tQ, Gl ) . 

Initial utttn.* Uj , iLT, iJbK- -©lyvatlon, gtot'raphicttl lati- 
tuaff, lorifc'ltua# ai U.a |<oinl 1, 

K«i»aitRj Blweit Xr<3/— Grtfciiwicu couralri»L«#u uf ti*i point n /bi> 

GP- twodi-M'.- ..’i lii ok (5, 5) i.atrix cf the tr'-nRf'uraat Ion froa 
ii/uiou a vooraii.ai**.. L(*. *...« ,cJa, rsyr.tuu., 

i, Uoe of the domain CCKIKfKj uno uaey tne conatanta fro". t..y 
biouk /GAK/~ leeu (,), Ho. li, fnole c‘,1). 

I 

7, Alg rlt m; Cooidtnatea of point i in lito 3r enwicu ay. tern 
of conralnate no dslermtned wltii tne formulas; 

^PO 9 

ypo *(<3te/4 •*-hp)cos(pp ^inKp, 

^PQ +Ap)*in ifp , 

A “ (cot*<pp + ( I 


nii ha^, /jI ; ha**-eomliBv.jor oxln and conf ■’icLer.l n comprea Ian af 
tliO tjarti.'U oiilpiiaio, 

Katrix of t£ie iraiuuuriaat* ;* from tao Grotnwioti 6,,ateni 

X gi 

of cuardlnato;! to lua puxiu. nyutom (aalornunou xn par, 1), v;uica 
hay li.o fuilo.‘. iuj form; 



»in (fpCcsKp 
CX)5 {pp cos Kp 
-sin Kp 


•in (pp sin A.p 

cos^p 

cos <pp sin Kp 

sincpp 

cos Kp 

0 



:+5 


o. Taxt*. 


SUIK0UT|N| POtHrlHP,ALTP0iLMP,Mp0tPt 

DIMINSION X^(SI.6^(3>St 

CO¥(iON/eAl/AC»ii 

ClteOluLT^I 

SlBSlNUiiT^I 

CbNaCOilALN^) 

SlNallNf AlNpi 

AiKlALK*ALK 

• N»f8«T*C|.»CL*AAK*ft*H> 

ABAf/AN 

6^1 ii n>-SL*CLN 

6A< |*tM*SL«tlM 
6AU,)).C|. 
saOi 1 >»ei*eiLN 


0A(}»|)««SLm 

SAO«tUei.N 

6AO«S)bO. 

VAf >}t( A*ALK*HP}*St. 
aba*ha 

MA(||«a«6A(2.1) 

)(A(})sa.6A(2«2> 

aituan 

(NO 


f iv /*’•». »Si it*f% i“v» k.rt I . 44 

|.8iW« * V-*^ t-ri f- 

t’.ipocoiit rlc form, dcitormi 
/iX.iitmth (U 3 — :mx, '»TK!) 


%iJ|S I »“a 

ii rr » V AS 

:r-t i os* 


#» '» »'C !• ^4*Vi P * t\ «•■ 

' JL i ii ‘ * V’ *^7 y- 

.•f ron^Lng, of locoii.jn, 


•t k4 ^ 

i,ix \f 


/• 




1 . iurro^.ie. .3ubprot:,i\'o,. :< for kno’AT. d?‘e* nwLcr coord Inatef? 
of Uu» AES y» 1 ') ona '’reenwio;. c. >rdln !»•« af tiio point I 
(x^ f,f jo t, ,..) 'letoruiiueo tno coord iiio too of abSh » jTj , ) In 

l^oint-tfji'OCQntric ayerom of cooi ainotesf IjAL, incllnedi rnnge 
(dtniaacQ Sr -u AES to . ini i); ^V, anglo ni Inc ti>-n--nnEle oetween 
rauiuB vector of AES, ana the fflarie, perpendicular to the axla T 

of litf p;iiot topuctutric syatem of coordin;»toc; A.i aaimutu—angle 

uetwe^n uiroctUm for i.!io aun (angle i , ) ana the pr.-jectlon of 

‘ , 

1 

the latilua vector of the AES on the |•lalu^ perpundicuiar tu the axle 
r.. (or-oaite the hand of toe cIoca). Tne componcmta of ti.c vnctur 

of velocity In the point aystem isf Ci.jordln-.tea may be dett** mined, 
uein thf* .'.■ubprngram^TJHN, intended for the raul tipi ic.' t Inn of too 
matrix of dimonsiun (N, K) by the N-dlmenslonal vector, 

A, 'itructure. The eubprograme GRfX, TURX, 

Sxternal oubprograme utilised! GCLTLN (RIO) 

3. Converalon to GHiX; 

C LL GRiX (XG, XGI', GF, Xi , DAL, AH, AH), 

Initial d-‘ta; Block XG 7 — Greenwich cooroin* ten of AES- block 

XGFj^ — Greenwich cnordln tea of point F; GI —twodimenslonel block 

(3,3)-“Tietrix of tho tronsformation fr m the Greenwich to the point 
syate of coordin»*ites. 


k) 



Rttuuiiaj block — coordln. tas of AES in topoctntric point 

ayatera of coordinates, DAL, AM, Akj, determined above, 

J+, Conversion to subprocrara TURN, 

CALL TURK fX, A, Y, K) , 

Initial dsi,a! X — block of dimension N; A— block of dimension 
(K, K) ; N— dimension, 


Results: Y— block of diraenHion N (Y=AX), 

6, Algorithm 




t 




^ a 

X ’- XpQ 



». 

«[GPfy] 

y - VpQ 

• 

» 


w •• 

• 

Z - ZpQ 

m ^ 


— 1 
a. 

i 



[>A'L = ((^ -Xpaf -f* iy-Vpof 

hM *= arctfi\Zp/(x]i^yp ) ) ; 

AZ « Arctg CypA/») , 



to the 
7 


-matrix of transformation 

i-. 'nt form; 

Text; 


from Groonwich system of coordinates 
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.SUBROUTINE CRPXIXU«x6P#6P#XP4pAi iAM»AZ} ^ 
dimension XGlS)»XCP($)i6P«'S#3l*xPl3l»XC$) 


OAisO. 
00 2 


X(J)sxc( J>»X6P< 
DHsCAl*X( J>*X< J> 
DALsSQPtIDAU 


CALL TU«NlX,CP»XPr3) 

xm*xp(U 

X(2)SXP(3) 

X(3)CXP(2) 

CALL GCLtLNtX* AM*AZ) 

RETURN 

ENO 


subroutine TURN(Xg#CPkXPiN) 
dimension Xe(N)f6P(NiN)»XP(W 
00 1 I=i»N 
XP<I>*0. 

00 1 Jsi*N 

I XP( I>sXpn><»GP< t»J>*X6U> 
return 


i 


end 



4.V. Transformation from the absolute system of coordinates to 
the solar-ecliptic ( 109-^-AH3BCL) , 


1. Purpose, For known coordinates of ASS (X, Y, :') in the 
ab'^olute system of coordinates one determines the coordinates of 
the ASS in the solar-ecliptic system of coordinates, determined by 
the following means: the center of the syatem 0 coincides with the 

center of the earth, axes CX__, OY.^ lie in the plane of the ecliptic 

the axis directed towards the sun, axis 0 perpendicular to 

the plane of the ecliptic and forming an acute angle woth the axis 
of rotation of the earth directed toward the north. 

B, Structure, Subprogram IBSKCL, 

Subprograms utiliaed; QGLTLN (BIO), Common block /BIECL/^, 


3. Conversion: CALL AB5ECL (XA, XAS, XSE, CLTS, ClNS). 

4. Initial data; Block XA, — absolute coordinates of IS”.; 

Block XAS^--directed cosines of radius vector nf the sun in the 

fibsolute system of coordinates. 

r. Hesults: Block XS37 — coordinates of AES ; CLTS, CLKS — 

geocentric latitude and longitude of AES in solar-ecliptic system of 
coord Lna ten. 


of Uie domain 


COMMfK. 


From the block C(rfCX /BIECL/- 


on© uses cos epsilon and sin epsil >n, computed in tht? subprogram 
SlIN(Cti-), whuro epsilon — angle of inclin- ti :ii of the plane of the 
earth's enuator to the plane of U’.e ecliptic. 

7. Ugorithui; 



8. Text: 
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8. Text: 


t' 


lttftflOUT|Nf AftSlCk(V»Xf*XfC*AIC«tiC) 
cevnOM /eifcL/ci^tSip 
DiMENStoM ViSIfItOl.XECOI 
M|asail?CXS(2)*XS(aUxSOUXtf}l) 
irixna.a.a 
NBCfa 
>Nl«SCa 
•OTO « 

NbXS(2)/XI 
Ml«XSt 3 )/Xl 
XECI IMe* 

DO 1 J«l *3 

xieiiivYief i)*v(4>«xt(ji 


XietSlayf 3l»M.V(Tr*M| 
lwaS|GN(i.*X|C2n 
XlCftlaxECCnal* 
XEec3)axEC(S)a|« 

CALL OClTLNCXECfAEC.eiCI 

aiTUXN 

END 


< “* » 


t.ateul tl 


4»T0. Mttiae.U af tUn- of the e.ary of <iES Into the e-irth'e shadow 
and exit. ir *m the .‘h.auow (110— .hi O'W) 

I, iar:u.;e. For known elements of the ■.i-'bit; a, e, i, 

Urae^a, ome^a, for rositiun in the orbit (u) at a ^iven moment in 
time t one dotermines the wumentu of time of entry of the I.lij Into 
the shade ’.V of the earth and exit from Ine suadow. 

c, structure; ."ubjsro^rem JilnhH.V, 

One uses external ,;abpic;/:rams I?uDT h (III), (GOF), 

Common blocks: /CGR/*^ , 7ci?ii)/ ^, /Cil/^, 

3. ConverRlnn: 0^1,1. dLADih'/ (DT, T, A, ?i, T2)» 

4 , Initial data; Block A^, contatninc elements of the orbit 
ami the ari:;unont of latitude: ^a, e, i, Gme^ia, omega, u; DT-~data 
in the form RJD, T — Moscow time, corresponding to tne argument of 
latitude u, 

!3. Results: T1 , TO — are, respectively, the moments of entry of 

tho ^ES into the shadow of the earth and exit from thesViadow (in units 
given by the scale factor EdEC). If AES does not get into the shadow- 
then T1^TS=v). 

6, Use of the domain COMMON: One uses the constants from the 

blocks CuiiMc N/CGR/^ , /CRE/, , /Cll/y (see (3), Nos. 5, 9, 

table 0.’, No. 1 table F.3). 

?. Algorithm: Moments of time (T1) of entry of aes into the 
shadow of the earth and (T2) of exit from the shadow are determined 
without computation of the effect of che penumbra j with the assump- 
tion tnat the earth has no compression and is not displaced orbit- 
al.! y. (7). These moments of time correspond to actual anomalies of 
the AES, satisfying the following ^‘elationship; 

cos ^ ^ ^ (1) 

where — radius vector of tho sun, r — radius vector of the 

AES, B.Q-- equatorial radius of the earth, pSi -angle between 

radius-vector of aes and radius vector of the sun. The relation- /6l 

ship (1) leads to the equation of the fourth degree with respect to 
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the cosines of the actual anomaly: 


^A,coa*v *AtCofv*A,eo^v*A,cMV*A^ , (2) 

>1,-6 (^)V- 2(^)’ ( 


+2a*-3’)(l-£*)-4S»4*. 



P-‘XlPj*YlPy*ZlPt, 

i ,^XqQs *^o^y • 

- HanpaBJUDDcme KocnHyca pawcapBeKTopa Cojmna, Rq , 

0 , 07 •» directed cosines^ of radius vector of the sun 

« COS CO cosQ “ sin to slnO cos i , 

Py - cos to sinQ + sin CO cos Q cos t , 

P^ - sin CO sin i , 

-sinco cosQ - cos to sinC^ cost, 

Q y* - sin CO sin Q + cos to cosQ cos i , 

0-» cosco sin t . 

a The solutions to equations (2), satisfying the conditions; 

beta cos v + xi sin v leas than 0, ^ _ 

S= (1 +ecos v)'^ + (p/a (beta cos v + xi sin v)^-(p/a =0. 

t? 6 

are assumed for further conaidei’atipn. 

The condition; 

j Smp\^cosv ^^sinv)(tcosv^f^sinv)-a\e(i4‘eco%v)5\nv >0 

determines the value of the actual anomaly v^ in the entry of theIS2 
into the shadov/ of the earth, S less than zero correspoinds to the 
value Vg, determining the exit of thelSZ from theshadow. 


53 


corresponding given 

thf«^ t°J„!!?.; tau, the tima of 

t - t-(C.-<*ln£,)/A, 

I rj» A •/ins' , ts£./2 .([’i-c)/fi.»c|'^teV2, 

constant times the mass of 

Moments of Lime (isl, c) are determined by the known 


Ti -r ♦(£^.e8ln£*)A , 

i^Ei/2 ■f(l-e)/(l+ej)*^t£p|/2 , 


8« Text: 


^UfROtfTf Nf •T»itnifll 

DtNIHSteN A|«|,8f«f*t«fl«fl4l»lftCt> 
oirtErftoM ifili 
eOMMAN/CM/ft . . 
eoHMfN/eAf/tt «f 

common /CAl/AttV 

Tl*|, 

Tl»l. . 

00 I J9%$j 
R«U.U«f4t 
OCKMOOilAUotll 
I OfK#l>«SlNtOU«t>l 

Ai.etittCAIt|/OOMA<|l 
V*(A|0|*0Cf II 
NtMfNCfl . 

wi«eoftv»*l,. 

trtMi.ot.i«i*t»i 

. VMfAN<Nf«fRmi«t. 


mf.iT.o.> 

V«0*0lt 

m«v*Atii*S|Hevi 

CiLl .SUNlOT,T*»S#*t,N»tttl 

trfocii,ft,oi 

OOTO S4 

Nta.xSlt}*0|4t#lt(Sl*0t0l 
N}«»xS(ll«0|4l«KSttl«0f 91 
HSe«r4NtlNt,i;i«0UI| 

VaM>*Al9l 

MAtSfMfVl 

Nt»eOtfVl*l,s 

tfCMt«eT*t*i«t9i 

. . V«ATANf«8«EK/Ni|*8« 

ir(v,iT*o«> 

. oaV^Oft 
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passage of 


tl.e earth, 
values t and 


4 


OEIGiNAL PAGE IS 
0P POOR QUAUTY 


C*Ll 

94 CONTtNUE 

AE*Eai*E(f| 

AQ«E(1>«E(?) 

EUM-PUl 
EI4>sPf 
DO 2 Jmifl 

PfJiaPf Tt*Pf.|*9l»AP*PIJ*4l 

PfilsPff I 
6*9 
$•0 

00 SJ«I»9 
0«XS(J»*PU »«9 
9 SiX$(J >*P| J*9>#S 

Wt>A(2)*A(2) 

PK«i(l)*a*.»*l)/4r 

PX«PK»PK 
W2«l ,/PK 
H1«W2*W2 
* AMBfB 
Afl«$«S 
W5saQ*AP 
W4ti ,-Ao 
P(«)s2. 

P<9)r4,«A(?) 

P(4)c2,«A(2»*A(2 
P(3 »s2.*P<4> 

PI2)*P<3»«A<2) 

P( I )»P(2>«Pf ?)/A 
Q( 1 )cAP*AQ 

Q(l)afi <*>•«<:> 

Q<2)s0 

!• »V3isq«T(i,-v( j) j 

WlsB«V(J> 

H5»i .«Af 2l*v(J » 

H4*W5*W9»PK 
00 36 111. 2 
H4«S*N3 
W4swi**’H 
V«H6«PP«W4*iu 
ircAescv»-i.!;-S)2^»30»90 
2'^ fPfw4)19*30,30 
19 V»PI 

IP(w2.CT.1.E»19) 

. V«4T4M(M3*EK/W2)«2 . 

IFCV.LT.O, » 

. VsV«Pl2 

VsT«( V.A(2)«S1P(V)-TP )*AL 

23 W2sPK *l^4*( S«V ( j )*B«W3 >-»«9«A f 2 >«h3 
IF (H2)24i29,23 

25 Tl*V 
goto ‘36 

24 T2»v 
30 W3*-i»3 
17 continue 

▼ PETtTN 
ENO 



0(3 >ft2»'#Tf3«Vr4.4,»AP»AQ 
a(4|«0 

Q(3)«W4«H4»w> 

DO 9 i/cl»3 

0(J)»Q(J)*Ni«P(^ »-»»2*PC J*3»*WJ- 
9 Q( J*2»*Q( J*2)*Wl*OIJ*3>-W2#PI J*3MW* 

CALL ROOT4(o,Vf IftI 
N«4 

K«l 

IF(|P(n)8»A»4 
6 KX3 

8 |FUC(2n9*9,ll 

11 N«2 

9 ir(K-N)l2«l2*7 

12 00 1? JtK.N 
trtAPSfufJIt.l tlAtlS*l7 
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4.1 it Calculat'on of the roots of algebraic equations of the 
fourth, third and second degree (I11-RC4T4) 

1, iurpose. Determinat iun of the roots of the algebraic 
equation 

A^x^‘ + + A^x^ +A^x +A^^ sO, 

of the fourth A^/ 0, third (A^-0) and second (A^^^u, A^-0) 

degreo with real coefficients. 

2, ijtructure: r>ubprugr.4n HoGT 4. 

3, Conversion; CALL HOOTq (A, Y, IR}, 

4, Initial data; Block A — coefficients of equation (1) in 
order; A^, A|, A^, A^, A^» 

3. I?esult8; Block roots of equation (1); 

Block IR^ — characteristics of the roots: 

0, if Y('l), Y(;!)--roal roots 

4» if Y(l), Y(,‘) real and imaginary part of a 
(1)= complex root (for complex-conjugate root ¥(2) 

is taken with the op-osite sign) 

1, if Aqj; A^r -^2"^ (equation of first degree) 

0. 1f Y(3). Y(/i)--,real rt,-ote 

4, if Y(3), Y(4)~“r9nl aiiij iinaginax'y pai’t of 
co;i!piex root (for complex conjugate root Y(4) ta 
IR (r) - kenwLtn optosite sign) 

3, if Aq=C (equati m of third degree) 

2, if Aq-A ^-0 (equnt'n of second degree) 

1, ii Ay=A^.-A^,^0 (equuwioii uf first dogr«e). 

Hejuarh... Fur IR (I) si Y(1) — root of equatinn -li Imj iir^v 
u'-gree. For IR (4) -3 Y(5) — rn„i root of an eouation of tiic thiru 

uogx'ee. 


algtxrt .1 Jim 

a ) i'j 

Dividing 


igoritl.. Fo, he solntiou of equation ( 1 ) one uses the 
, described in (7) (p, 446); 

quation of 4 degree is solved by the method of Descartes: 
equation (1) by Aq (if Aq- 0, then we proceed to part delta). 


we obtain 




* 0 , 




V»e transform thisecuation, excluding from it the term of the third 
degree: 

+Qje +,# « ^ 

^ +B3 , 

h ~ -B| /4 f j 9 ^ ^ h , 


po 



For '*=0 m obt?iin the biquadratic e-iuatiunj for the solution of 
which we uoe the formulas of p rt v. For J 0 we consider the 
cubic resolvent 

+C,t mO , i 


c,»ae, c, .-o%| 

I 

We transform equation (3) into the form: 

z ® ■♦■A -0* 

O /3 «C,/3 - 3 ‘ ^ /2 « ♦«•- e, 3/2 4 cv/a , 

« * Cj/S , 

* _ a 


» «t43 , 
* 
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(3) 


(it) 


ir lJi;lta is greater than 0, than the solution of eouaiion (4) 
is found F: ■*“ rnrwuia of Kuolan: 

Z.,=(-6/2 +(A)''*)"'* *(6/i-(A)*'*y'*\ 


I • Xs - KomneRcme RopBi. 

z^f E^^-'-compiBA ruuvto. 


If Delta =0, then the ro ts of equation (4) are determined by the 
formula ; 


1 A sO, TO RopHB ypasBeBJBA (4) onpeAejnen no c 

U Z.mZ. 


If Delta is less than‘0, then we use the representation of the 
roots as resolvents , in trigonometric form, 


•EoCOS<(>/i , 


B £oCOs(4>/5 4 120*) , * £«cos((^/3 4 240^ 


i 


We select a critical root of eouation (3); 

N. 

where z^-s— real root of equation (3). 

Knowing R', it is possible to break down eouation (2) with the 
multipliers; 


(X 

I' 





Usinc tiiereliitionshipe of division V, w© obtain the solutions 

(i=1,4) of two quadratic snuations (i), and consoouently — th& roots 
of squ-*t i'm (1): y. =x^+h. 

6. Cubic equation we write in the form 

y^ + A.y^ + A^y •U^^=0. 

Dividing the left side by (if A^sOf we go on to division b), 

we obtain ^3 

Tnio equation coincides with equation (5), the solution of 
which is described in division a, 

Humark For Delta greater than zero the complex roots y^ , y^., 

of equation (u) are obtained by the solution of a quadratic equatitn 

y • + (c, + *, -*)y + e, + (c, 4 -a)( * , -j) , 

whore a^— roal root of equation (q), ssc^/3. 

b) Dividing the left side of the quadratic equation 

4*« 

o^y*" +A^y +A^iO by A^ (if A|j-0, then we go on to 
division c) , we obtain: 
i +P,y +B,^.0. 

The roots o<‘ thisenuatlon 


y,., . 


c) Fnuati »n of the iirst aegree A^y +A^^-0< 


Solution y~ “A, /A* 

'i J 
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7, Text; 

OtMCNStoR *( 9 >»VU)* 

DATA P|,RI0|,9|*/9.iAltftA9,l.fTi.tfRSI»i*tM|iln/ 

CATA |N/1»9*2/ 

K«i 

IRtlUO 
« tRttMl 

Iff A(l))l*ll9*i 
i K*A 

DO 2 J*|*A 

i »l«l»BAOfl)/A(l» 

DO 3 

3 V(2MVC2>«H*D(J) 

a»A. 

*»•>- 

DO A tlBlil 
Q«ReH«d( J)*M 

4 MeM*l. 

pe(2.*H*DU)»*3.«H*d(2) 

VU)«P 

N*l 

lMADS(Q)-l.e-lAU03»i03*lOA 
109 00 7 J«l»Ni2 
|BlNU> 

0BV(J)«V(«t>-4,*V(il41) 

J)*-V< J W2, 
lr(0>9*6*6 
0 0 «SqRT(0)/2. 

.V(J«U*V(J>*0 
V< J)*V< J>^0 
tR< t )*0 
OOTO 7 

9 Vi J«l)sSaRT(*0)/2. 

7 COiiTlNUl 

60T0( 104* 10»« 167 > *K 
104 lr< 1R( I ) llOillt lO 
11 ViA)B*V(l) 

V(2)»-V(2) 

V(9)«0. 

v(n>o» 

•N»? ■ 

K«9 

OOTO 103 
10 Ml«v(l>*V(l) 

W«Wl*V{2>*V(2) 

RIBSQRT 

W»SqRT(W*WI) 

R1»SQRTIR1>*2. 

V(1)*V(2)/R1 
V(2)««i/Rl 
V(3)SH»yill 
ya)*v(i»*H 

W(4>*V(z) 59 

lR(2)s4 

• 1« •». ». 


PAGP ic 


* 


6B 


60 


100 ir U<2> )19rlfll»l> 

10 Do lO 

16 6( J»»A( j*l> fk^7.^ 
IR12)»3 
OOTO lot 
106 »(2)>VUM?. 

6(3>rv(i>*V(n-«**V(21 

eiA 

iC# 5*Bl2 >/J. 


V^kKaS 

ARRftOi/S.-m 
6R*S*6I ?)/* .-5»W-8C A»/2. 
prAp>*AR.AR*6ft«6R 
IFCD) 17,16,10 
i? F«k2.«S0RTI-AIU 
11F({}R>20p21 ,20 


21 HlsplOZ 
00T02J 

20 M»:saRT(-D» 

|F (6R>5l»9fl»9l 

50 Wi»pID2 
SOTO 23 

51 fil*ATAM(Wl/90> 
1FC»R)22»23,23 


II srirsi.-j 

W<Jj«W.cOSfWl/S,l 
2A N1«w1*Pi 2 
' GOTO 110 
19 WsSeRT(O) 

W<1)*BR*W 

VI2>*B0-W 

Do 29jBi,2 
MxAoS(V(<n ) 
wi*aloocwi 
WI*EXP(wI/9. } 

29 VI J )*SlG»<<Wl, V( JU 
VI3| »VI1 >*VC2)-S 
y(i)ce, 

•V<2)*0. 

X»IR(2) 

hl*V<3» 

GOTO ( 111* 209,209) ,K 
18 pftoesieA) 

‘ Wl»AL0C|Pn 

W1S6XPIW1/3, ) 

V ( 1 )«SIgNIw 1,BR) 
V(2)B-y(D 


107 
; >l 

i, ’» 

;!*o5 


y ( 1 ) b2«v( 1 ) 

Vt3)sV(2) 

110 0C26Js1,3 
26 V(J)*V(j)-S 

GOTO ( 106, 109,209) ,K 
106 risamax 1 <V< i),y(2),y(3)) 

111 V( I )sSGrTCR1) , 

Vl3)»-V(l) 


icl 

36 


WsPfRl 


WlsQ/V(l ) 

y(2)MM.Wl)/2, 

V(4)clM4MD/2. 


102 

37 


K«S 

GOTO 103 
DO 30 Jsl,M,2 
ISINI J) 

VC J )»V| j )*H 

IM IRC )t 5)30,31.30 

V ( J*l lsy< J*l >^M 

CON?lfiUg 

CfOTO 109 

yCl)s6(2l*V<3) 

V(2)SV( 1>»V<.3»*6<3) 

N*l 

Ks2 

GOTO 103 

IFCa<3) >36, 102,36 

yCl)=A{A>/A(3) 

V(2)sA(5)/A(3) 

1R(2)«2 

Nsl 

GOTO 103 

lr(A<A))37,lo9»37 
V(l)>-Af 9 )/a{ A) 

GOTO 109 
END 




% 


UZl F nUP4'm aR/O-iri ISY N^.F.E 


NAME 

— TJiiJKX 

' f> r* L* 

NAME 

Y *»tH1H** 
. *1 

! *GS 

A&SECL 

./ 09 

58 


LIMES 

/ M 

52 

AmMSP’ 

QOI 

24 


LOQMOD 

F xn 

18 ‘ 

ADWr 

002 

31 


MEWMAO 

I 04 

•92 . 

MNIAT 

. / 05 

52 


MODE 

H 03 

41 

APSIS 

H 04 

43 


OIMLAT 

I 05 

5? 

ASfN 

/ 04 

52 


PARAB 

H 02 

. 36 

AULOm 

I 03 

50 


POIMT 

/ 07 

55 

BL 

/ 04 

52 


PTERM 

I 06 

54 

CARMEL 

I 04 

52 


REACH 

H 01 

33 

EXTRM 

H 02 

36 


RIQHTA 

I 01 

46 

FAORAV 

F 03 

II 


RODENS 

F or 

18 • 

EATM . 

F 02 

10 


ROOT4 

I II 

64 

ATCllOO 

F 06 

17 


SHADAW 

/ 10 

60 

PORSS 

F 04 

■ 14 


STAR 

/ 04 

52 

FLIGHT 

F 05 

15 


TAUR 

/ 03 

50 

FNQRAV 

F 01 

8 


TOEOM 

I 02 

48 

FORCE 

F(T7 

18 


TLOC 

I 01 

46 

OMLL 

I u2 

48 


TURN 

I 08 

• 57 

ORPX 

1 08 

57 


ZDSUNC 

/ 01 

46 

INTEO 

/ 04 

52 


ZbSUNL 

/ 01 

46 

fWAA 

‘ / 04 

52 
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